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Abstract: We propose a method for matching the next-to-leading order (NLO) 
Oh. calculation of a given QCD process with a parton shower Monte Carlo (MC) simu- 

^ , lation. The method has the following features: fully exclusive events are generated, 

with hadronization according to the MC model; total exclusive rates are accurate to 
^ . NLO; NLO results for distributions are recovered upon expansion in as] hard emis- 

sions are treated as in NLO computations while soft /coUinear emissions are handled 
by the MC simulation, with the same logarithmic accuracy as the MC; and match- 
ing between the hard- and soft /collinear-emission regions is smooth. A fraction of 
events with negative weight is generated, but unweighting remains possible with rea- 
sonable efficiency. The method is clarified using a simple toy model, and illustrated 
by application to the hadropro duct ion of W^W^ pairs. 
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1. Introduction 

The reliable prediction of cross sections and final-state distributions for QCD pro- 
cesses is important not only as a test of QCD but also for the design of coUider 
experiments and new particle searches. All systematic approaches to this problem 
are based on perturbation theory, usually truncated at next-to-leading order (NLO) , 
or in a few cases next-next-to-leading order (NNLO). Such calculations yield the 
best available results for sufficiently inclusive observables, in kinematic regions where 
higher-order corrections arc not enhanced. However, in many cases a more exclusive 
description of final states and/or a wider kinematic coverage are needed. 

For the description of exclusive hadronic final states, perturbative calculations 
have to be combined with a model for the conversion of partonic final states into 
hadrons (hadronization) . Existing hadronization models are in remarkably good 
agreement with a wide range of data, after tuning of model parameters [1, 2]. How- 
ever, these models operate on partonic states with high multiplicity and low relative 
transverse momenta, which are obtained from a parton shower [3, 4] or dipole cas- 
cade [5] approximation to QCD dynamics and not from fixed-order calculations. 

To extend the coverage of perturbative calculations into regions where higher- 
order contributions are enhanced, one needs to identify and resum the enhanced terms 
to all orders. The resummed expression can then be combined with the fixed-order 
result, after subtracting resummed terms of the same order to avoid double-counting. 
This has been done explicitly for many important observables and processes, but a 
general procedure has not yet been developed. In principle, the parton shower pro- 
vides such a procedure, if the problem of double- counting can be overcome. Combin- 
ing NLO and parton-shower computations also offers the possibility of interfacing to 
hadronization models, thus overcoming the main deficiencies of existing calculations. 

Methods for combining fixed-order and parton-shower or dipole-cascade compu- 
tations in QCD have been discussed in refs. [6]- [21]. Many of those papers concern 
the rather different problem of combining leading-order (tree level) matrix elements 
for hard multi-jet production with parton showers or dipole cascades for jet frag- 
mentation. The only fully-implemented schemes for matching NLO computations 
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and parton showers that have been used to produce phenomenological results [12]- 
[15] have been based on the so-called phase-space slicing method for regularising the 
infrared divergences encountered at NLO. The slicing method is simple but has the 
disadvantage of introducing a slicing parameter, which should in principle be sent to 
zero at the end of the calculation but in practice cannot be chosen too small. This 
can introduce uncontrolled errors and/or large negative event weights. The slicing 
method has therefore tended to be replaced where possible in NLO calculations by 
subtraction methods, in which no approximations need to be made and event weights 
can be controlled by improving the choice of subtraction terms. 

The subtraction method has been considered in refs. [11, 18, 19, 20]. The main 
concern of those papers was the correct matching of NLO matrix element calculations 
to the scale and scheme for next-to-leading logarithmic (NLL) parton showers. Up to 
now, no satisfactory algorithms for NLL showering have been developed. When such 
algorithms are available, the methods proposed there will make it possible for event 
generators to generate the NLO evolution of parton distribution functions (PDFs) 
and hard process cross sections, without the explicit input of NLO PDFs. Existing 
parton shower generators, however, correspond only to LL evolution, and therefore 
the matching criteria of refs. [11, 18, 19, 20] operate beyond the level of precision 
that can be achieved at present. 

In the present paper, we present a method for matching NLO matrix elements 
and existing parton shower generators based on the subtraction method for NLO 
calculations. First, in sect. 2, we outline the general objectives and limitations of 
our approach. The method is then illustrated in sect. 3 in the context of a simplified 
toy model. The case of QCD is dealt with in sects. 4 and 5. In sect. 4 results are 
presented for the process of W"'"W~ production in hadronic collisions. In sect. 5 
we discuss briefly the application of the method to other, more complex, processes. 
Concluding comments are given in sect. 6. All the technical details are collected in 
the appendices, where we also comment on the matching between NLO computations 
and parton shower approaches in the framework of the slicing method. 

2. General framework 

2.1 Objectives 

It is undeniable that parton shower Monte Carlos (MC's) provide more general and 
flexible tools compared to NLO computations. By running an MC we get events, 
which we believe to be a faithful description of the real events detected by real 
experiments; this is one of the reasons why the MC's are so popular with experimen- 
talists. On the other hand, NLO computations emphasize the role of (infrared-safe) 
observables; the presence of negative contributions, and the necessity of considering 
inclusive quantities in order to get rid of infinities, make it very problematic to even 
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talk about events. Still, NLO computations have their virtues; they can handle hard 
emissions, and they can estimate total rates with better accuracy than MC's. 

In this paper, we aim to construct an MC that works like a regular MC, but on top 
of that knows how to treat hard emissions, and can compute rates to NLO accuracy. 
For brevity, we call this tool an MC@NLO. Clearly, there is a lot of freedom in 
constructing an MC@NLO. We shall follow a minimal approach. Namely, we want to 
keep the original MC as unmodified as possible; all major changes have to be carried 
out on the NLO codes. Amongst other things, this implies that our MC@NLO 
will not improve the logarithmic accuracy of the original MC, in those regions of 
the phase space where resummation is needed. As far as the NLO calculations are 
concerned, we shall not perform any approximations, neither at the level of matrix 
elements, nor at the level of phase space; thus, we shall adopt the subtraction method 
in order to cancel the infrared divergences that arise in the intermediate steps of the 
computation. 

To be more definite, we require MC@NLO to have the following characteristics: 

• The output is a set of events, which are fully exclusive. 

• Total rates are accurate to NLO. 

• NLO results for all observables are recovered upon expansion of MC@NLO 
results in as- 

• Hard emissions are treated as in NLO computations. 

• Soft/coUinear emissions are treated as in MC. 

• The matching between hard- and soft-emission regions is smooth. 

• MC hadronization models are adopted. 

In sect. 4, we shall construct an MC@NLO fulfilling these conditions. However, 
we find that the technical comphcations of QCD tend to hide the basic ideas upon 
which our approach is based. Thus, we prefer to start from an oversimplified (and 
unreahstic) case, where the structure of the MC@NLO is apparent; this is done in 
sect. 3. The reason for doing this is that we shall see that the approach that leads 
to the MC@NLO in a simple case can actually be adopted in the QCD case as well, 
which will only be technically more involved. 
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2.2 Negative weights 

We do not require the MC@NLO to be positive definite; events may have negative 
weights. However, it should be pointed out that these negative-weight events are of 
a completely different nature from the negative contributions that appear in NLO 
computations. In particular, the distributions of the positive- and negative-weight 
events are separately finite in our MC@NLO (the same is not true in an NLO compu- 
tation). This implies the possibility of unweighted event generation,^ as is customary 
in ordinary MC's. The situation is actually completely analogous to that occurring 
in polarized collisions, where physical observables are obtained as the difference of 
two positive-definite and finite quantities. Furthermore, we shall see that in our case 
the number of negative-weight events is reasonably small, and thus the number of 
MC@NLO events necessary to get a smooth distribution is comparable to that in 
ordinary MC's. 

2.3 Double counting 

When the MC generates events with real parton emission, it generates kinematical 
configurations that are also taken into account by the NLO computation: the possi- 
bility of having the same kinematical configuration from the MC and from the NLO 
may lead to double counting. We shall adopt the following quantitative definition of 
double counting: 

An MC@NLO is affected by double counting if its prediction for any observable, at the 
first order beyond the Born approximation^ in the expansion in the coupling constant, 
is not equal to the NLO prediction. 

According to this definition, double counting may correspond to either an excess 
or a deficit in the prediction, at any point in phase space. This includes contributions 
from real emission and virtual corrections. Generally speaking, the MC fills the phase 
space efficiently only in the soft and collinear regions. In these regions, the leading 
behaviour of the kinematical distribution of one-particle emission is given by the 
Born cross section, times a kernel that describes the soft or collinear emission. This 
holds both for the MC and for the NLO. Thus, in the soft and coUinear regions of 
the one-particle-emission phase space, the MC and NLO results coincide when only 
the leading terms are considered. Elsewhere, the MC is not reliable and the NLO 
provides the best estimate. The problem is to merge these two descriptions. 

^In this paper, we adopt a slightly different notation from that commonly used by experimen- 
talists. Namely, we call unweighted events any set of events whose weights are all identical, up to 
a sign. In other words, the weights of our unweighted events are of the form ±w, where w is a 
constant. 

^By the Born approximation we mean the lowest order at which the hard process contributes to 
any observable. 
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3. Toy model studies 



In this section we study a toy model that allows a simple discussion of the key 
features of the NLO computation, of the MC approach, and of the matching between 
the two. We assume that a system can radiate massless particles (which we call 
photons), whose energy we denote by x, with < x < Xs < 1, Xs being the energy 
of the system before the radiation. After the radiation, the energy of the system is 
The system can undergo several further emissions; on the other hand, 
one photon cannot split further. 

In a perturbative computation, the Born term corresponds to no emissions. The 
first non-trivial order in perturbation theory gets contribution from those diagrams 
with one and only one emission, with either virtual or real photons. We write the 
corresponding contributions to the cross section as follows: 

i;^)^-BH.). (3.1) 

'4) =«(^ + vVw. (3.2) 



dx J ^ \2e 

for the Born, virtual, and real contributions respectively; a is the coupling constant, 
possibly times a colour factor; B and V are constant with respect to x, and 

\imR(x)^B. (3.4) 

a;— >0 

The constant B appears in eqs. (3.2) and (3.4) since we expect the residue of the 
leading singularity of the virtual and real contributions to be given by the Born 
term, times a suitable kernel. We take this kernel equal to 1, since this simphfies 
the computations, and it is not restrictive. Finally, e is the parameter entering 
dimensional regularization in 4 — 2e dimensions. 

The task of predicting an infrared-safe observable O to NLO accuracy amounts 
to computing the following quantity 



(O) =lim / dxx-^'0(x) 



da\ / da\ /da 
dx)^ \dx)^ \dxj^^ 



(3.5) 



where 0{x) is the observable as a function of x, possibly times a set of functions 
defining a histogram bin. By "infrared-safe" here we simply mean that the integral 
in eq. (3.5) exists. The main technical problem in eq. (3.5) is due to the presence of 
the regularising parameter e; in order to have an efficient numerical procedure, it is 
mandatory to extract the pole in e from the real contribution, in this way cancelling 
analytically the pole explicitly present in the virtual contribution. 
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3.1 NLO: slicing versus subtraction 

Within the context of this toy model, the usual methods of cancelling the pole in e 
are as follows. In the slicing method, a small parameter S is introduced into the real 
contribution in the following way: 



{0)^= [ dxx-^'0{x) 
Jo 



i 



do 



+ / dxx^'^'^Oix 



do 
dx 



(3.6) 



In the first term on the r.h.s. of this equation we expand 0{x) and R{x) in Taylor 
series around 0, and keep only the first term; the smaller 5, the better the approx- 
imation. On the other hand, the second term in eq. (3.6) does not contain any 
singularity, and we can just set e = 0. We obtain 

{0)^ = aBO{0) J\x^ + J\x0{x)(^^^ + 0{6) (3.7) 
= a (^-^ + log 5^ BO{0) + aj\x ^^""^^^"^^ + o{5, e). (3.8) 

Using this result in eq. (3.5), we get the NLO prediction for O as given in the shcing 
method: 



(O),,^ = 5O(0) + a 



(51og5 + 1/)O(0)+ / dx 
Js 



0{x)R{x) 



X 



+ 0{6). (3.9) 



The terms 0{5) cannot be computed, and are neglected; in general, one must check 
that these neglected terms are numerically small, by plotting {O)^^^^^ versus 5 in a 
suitable range in 5. 

In the subtraction method, no approximation is performed. One writes the real 
contribution as follows: 



(O)j^ = a5O(0) / dx 



'2e 



X 







X 



l+2e 



(3.10) 



The second term on the r.h.s. does not contain singularities, and we can set e = 0: 

0{x)R{x) - BO{0) 



{O),^-a^O{0) + aj\x- 



X 



Therefore, the NLO prediction as given in the subtraction method is: 



(O)^^, = BO{0) + a 



VO{0)+ [ 
Jo 



dx 



0{x)R{x) - BO{0) 



X 



(3.11) 



(3.12) 



We rewrite this in a slightly different form, which will be more suited to discussion 
of the matching with the MC: 

•1 



{0\ 



dx 



Oix)^ + OiO) (B + aV-^) 
X \ X J 



(3.13) 
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3.2 Toy Monte Carlo 



In a treatment based on Monte Carlo methods, the system can undergo an arbitrary 
number of emissions (branchings) , with probability controlled by the Sudakov form 
factor, defined for our toy model as follows: 



A(xi, X2) = exp 



X2 

—a I dz 



(3.14) 



where Q{z) is a monotonic function with the following general properties: 

< Q{z) < 1, limQ(^) = 1, limQ(^) = 0. (3.15) 

z— »0 2— >1 

Note that the MC is also well defined without the last condition in eq. (3.15). We 
discuss later the specific functional form assumed for Q. li Xg is the energy of the 
system before the first branching occurs, then A{x, Xs) is the probability that no 
photon be emitted with energy z such that x < z < Xg- At fixed x, this probability 
tends to 1 when a —>■ 0; that is, no emission is possible in the zero-coupling limit. 
On the other hand, when a — > 00, we have A{x,Xs) — > 0, which means that it 
is impossible that the system does not emit. The toy MC code we have written 
implements ordered emissions in the photon energies^ and proceeds through the 
following steps: 



0) Define xq as the lower bound on the photon energy after a branching; thus xq 
plays the same role as the cutoffs used to define the Sudakov form factor in 
QCD MC's. 

1) Define Xm to be the maximum energy available to the photon; before the first 
branching, we have Xm = Xg, where Xg is the energy of the system. 

2) Pick a random number Ri e [0, 1] , and solve for z the equation 

)a 
= R^ =^ z^XmRi''. (3.16) 



exp 



dz' 



3) Pick a random number R2 G [0,1]. If Q{z) < R2, let Xm = z, reject the 
branching, and return to step 2). 

4) li z < Xq, reject the branching and exit, li z > Xq, a photon with energy z has 
been emitted. We let Xm — z, and iterate the procedure going back to 2). 

^This may violate energy-momentum conservation, as is usual in a leading logarithm implemen- 
tation. 
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This procedure ensures that the photon emissions have the distribution prescribed 
by the Sudakov form factor in eq. (3.14). As can be seen from that equation, the 
probabihty of having soft emission (i.e., with z < Xq), and thus of terminating the 
shower, is A{xo,Xs), which tends to one when a — > (weak couphng), and when 
Xs — > Xq (no energy left for non-soft branching). In this toy model, the hardest 
emission is always the first one. 

3.3 Matching NLO and MC 

In order to introduce some useful notation, we start from the LO perturbative result, 
which does not pose any problems, since it is in fact what is implemented in ordinary 
MC's. At the LO, the system has not yet lost any energy to perturbative photon 
emission. Thus, the maximum energy available to the photon in the first branching 
is 1. For each event, we set Xm = 1, and proceed as described in sect. 3.2. The total 
exclusive rate will have to be equal to B, according to eq. (3.1); thus, if we generate 
N events, each event has a weight equal to B/N. 

We can formally read this procedure from eq. (3.13), by setting a — and 
performing the formal substitution 

0(x) I^aiO,xUx)), (3.17) 

where I^ic{0,Xj,i) stands for interface-to-MC, and indicates symbolically the distri- 
bution in the observable O as obtained by running the MC starting from a given 
Xm, and giving each event the weight 1/iV; its precise definition is given in app. B.l. 
For a; = we need = 1; we do not need to specify Xm{x) for other values of x at 
present. Interfacing LO with MC, we thus obtain 

^) =S7mc(0,1). (3.18) 

/ MCOLO 

This equation reminds us of the total rate (B), and the maximum energy available 
for photons at the first branching {xm = !)• 

3.3.1 Naive subtraction 

The first attempt at extending eq. (3.18) to NLO amounts simply to substituting 
eq. (3.17) into eq. (3.13). We get 



da 



dx 







^ ,^ , ,,aR(x) ^ ,^ ^, aB 

X \ X 



(3.19) 



dO 

This equation suggests the following procedure: 

• Pick at random < x < 1. 

• Generate an MC event with Xm{x) as maximum energy available to the photon 
in the first branching; attach to this event the weight w^v — aR{x)/x. 
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• Generate another MC event (a "counter-event") with Xm — 1; attach to this 
event the weight Wct = B + aV — aB/x. 

• Repeat the first three steps N times, and normahze with 1/N. 

Unfortunately, this procedure is bound to fail, since the weights Wev and Wct diverge 
as a; — > 0. A small cutoff p can then be introduced; as a; — > p, we would have 
Wbv — ciB/ p and Wct — —aBj p. Both these weights are very large in absolute value; 
this means that the eventual unweighting procedure, performed by the MC, will be 
highly inefficient. Furthermore, as we shall see in sect. 3.4, eq. (3.19) has problems 
of double counting, defined according to the criterion in sect. 2.3. 

3.3.2 Modified subtraction 

To avoid the problems outlined above, we define the modified subtraction method by 

a[R{x) - BQ{x)] 



dOj^,^^ Jo 



dx 



{0,Xm{x))- 



+I^c{0,l)iB + aV + 



X 

aB[Qi,r) - 1] 



X 



(3.20) 



This amounts to saying that, in eq. (3.19), we subtract and add the quantity 

,^ ^aBQ(x) 
/Mc(0,a;M)^^, (3.21) 

X 

using Xm — Xm{x) in the first and Xm = 1 in the second term introduced in this 
way. We point out that Q is the same function that appears in the definition of 
the Sudakov form factor, eq. (3.14). By virtue of the second condition in eq. (3.15), 
the weights Wev and Wct attached to events and counter-events are now separately 
convergent as a; ^ 0. 

The two terms involving Q{x), that we inserted in eq. (3.20), are not identical; 
therefore, the procedure which leads from eq. (3.19) to eq. (3.20) is not a subtraction 
in the usual sense of an NLO computation. Still, these two terms do not contribute to 
the observable O at C(a), because they are compensated by analogous terms due to 
the parton shower BluciO, 1); this will be shown with explicit examples in sect. 3.4, 
and in general in app. B.3. We point out that the function aQ{x)/x that appears in 
eq. (3.21) docs not depend upon the MC cutoff Xq. It follows that the dependence 
upon xq of the MC@NLO, eq. (3.20), is identical to that of the ordinary MC, being 
entirely due to /mo and is therefore power- like. 

It should be clear from the description of the procedure given above that the de- 
pendence upon O in eqs. (3.19) and (3.20) is purely formal, serving only to allow us 
to write such equations in a compact form. We stress that, in the actual implemen- 
tation of all the MC@NLO's described in this paper, the computation of observables 



10 



proceeds exactly as in standard MC's: namely, events are generated, showered, and 
hadronizcd without reference to any specific observable. 

Equation (3.20) is our master equation for the toy model. We shall show in 
what follows that the MC@NLO defined in this way meets the requirements given 
in sect. 2.1. In sect. 4, we shall use eq. (3.20) to construct a QCD MC@NLO by 
analogy with the toy model. 

3.4 Observables 

3.4.1 Exclusive observable 

As an example of an "exclusive" observable in the toy model, we shall consider the 
quantity 

y = max{xi,...,Xn), (3.22) 

where n denotes the number of emissions in one given event (in the case of no 
emission, the event is not used to fill the histograms). This observable is the energy 
of the hardest photon of the event. It is particularly convenient because in the MC 
approach it coincides with the energy of the first photon emitted in the shower; 
when MC@NLO is considered, y is the larger of the energy of the photon emitted 
perturbatively, and that of the first photon in the shower (if any). This fact allows 
us to obtain analytical results for MC, NLO and MC@NLO. 

The NLO prediction for this observable is (here and in the following, we consider 
only y > xq in order to simplify the analytical results) 

^) =«^. (3.23, 
while the MC result, normalized to the Born cross section, is 

^) =„B«MAte,l), (3.24) 

dy J MC y 

which is just the differential cross section at x — y, computed using eq. (3.14), times 
the probability that there is no harder emission at x > y. 

In order to study the MC@NLO results, we have to consider the quantities 
luciy, 1) and I^icXy^^Mix)) that will be substituted into eq. (3.19) or (3.20). In the 
case of lyiciy-i 1)) the hardest photon comes from the MC and we have 

/Mc(l/,l) = a^A(y,l). (3.25) 

y 

In the case of hiciy, x^^^x)), the hardest photon may come from the NLO contribution 
or the MC, and so we find 



(y,XM(x)) = A^min{x,XM(x)},XM(x)^5(y - x) 



+a^A(y, XM{x))Q{y - x)Q{xm{x) - y) . (3.26) 
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As a shorthand notation, in the following we shall denote the contributions to da/dy 
from hiciy, 1) and I^^ciy , x^i{x)) by da^/dy and da^/dy respectively, since they cor- 
respond to a Standard MC evolution (i.e., with no prior NLO real emission), and to 
a Hard MC evolution (whose initial condition includes one NLO real emission) . This 
notation will be extensively used from sect. 4.1 on. In the case of the naive subtrac- 
tion method, eq. (3.19), these contributions are separately divergent, but their sum 
is not. We find 



da^ 
dy 



Q{y) /K( 1^ 
a A y, 1 

y 



^ dr 

B + aV-aB I — 

^ . 



(3.27) 



The results for da^/dy depend somewhat on the form of the function Xm{x) intro- 
duced in eq. (3.17). We define x,, to be the solution of 

and we assume this solution to be unique. Wc further assume x^i[x) to be a mono- 
tonically decreasing function; thus, x < Xi^{x) for x < Xe- Then: 
o for y < Xe : 



dai 
dy 

o for y > : 



a A(y,XM(y)) +a 



y dx 

— R{x)A{y,XM{x)) 



X 



(3.28) 



dfrn 
dy 



Ri.y) 2 



y 



Expanding in a, we find 

da' 



da\ 

dy), 



dy 



dy 



o(v) r^i^y'^ dx 

^ / -R{x)A{y,x^{x)) 

y Jo X 



-[R{y) + BQ{y)]+0{a^) 



(3.29) 



(3.30) 



Therefore the naive subtraction method suffers from double counting, according to 
the definition in sect. 2.3, especially at small where Q{y) 1. 

For the modified subtraction method, eq. (3.20), we find instead that both con- 
tributions arc finite: 



da§ 
dy 



= a A(y, 1) 

y 



B + aV + aB 



f 

Jo 



dx 



Qix) 1 



X 



(3.31) 



while 

o for y < Xe 



dau 
dy 



-[R{y)-BQ{y)]A{y,x^{y)) 

y fjq* 

+ a^^l —[R{x)-BQ{x)]A{y,x^{x)) 



y 



X 



(3.32) 
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o for y > : 

rdp\ ^ a^^^^^ _ ^^^^^^ ^ ^^QM dx^j^^^^ _ j^Q^^^^^^y^ _ 

V y / msub y y J ^ 

(3.33) 

Expanding in a, we now get 

^) =^+0(a^), (3.34) 

^ / msub ^ 

in agreement with the NLO prediction for all values of y. Furthermore, for y — > 
we have the expected rcsummation of leading logarithms, 

da\ ( aB , ^ 



^ . ( — + C»(a^) 1 A(y,l) + PST fory^O, (3.35) 

y J msub \ y / 

where PST stands for power-suppressed terms, i.e., terms that arc not singular for 
y — s> 0. We stress that the behaviour of da/dy in eq. (3.35) is dictated by eq. (3.31), 
i.e., the contribution of EI events is power suppressed for y — > 0, since {R{y) — 
BQ{y))/y = 0{1) in this limit. We also notice that the functional form in eq. (3.35) 
is the same as that in eq. (3.24). Thus, the MC@NLO resums large logarithms of y 
in the same way as does the MC. A general derivation of this result for a wide class 
of observables is given in app. B.4. 

3.4.2 Inclusive observable 

As an "inclusive" observable, we shall consider the fully inclusive distribution of the 
photon energies; thus, each photon emitted, either at the NLO level or by the MC 
shower, will contribute to the physical observable, which we shall denote by z. The 
NLO and MC results (for z > Xq) are: 

d(T\ R{z) 

— = a , 3.36) 

pi ^aB^^. (3.37) 

Eq. (3.36) coincides with eq. (3.23), which is a trivial consequence of the fact that 
only one emission occurs at NLO. On the other hand, the MC result in eq. (3.37) 
differs from the exclusive case, eq. (3.24), since all emissions now contribute, not 
only the first. As a consequence, eq. (3.37) has a logarithmically-divergent integral 

at 2; ^ 0, whereas this does not happen in eq. (3.24). 

In order to obtain the predictions of the MC@NLO, we need to find the analogues 
of eqs. (3.25) and (3.26), in the case of the variable z. We obtain 

7Mc(^,l)=a^, (3.38) 

z 

Imc{z, Xm{x)) ^ 5{z - x) + a ^^^^ Q{xm{x) - z). (3.39) 

z 



13 



The term with the (5-function in eq. (3.39) accounts for real emission at the NLO 
level, which always contributes to the inclusive spectrum when it occurs; the other 
term is the contribution of the MC shower. 

To compute the z distribution for the naive subtraction method, we use eq. (3.38) 
and eq. (3.39) in eq. (3.19), to obtain 



/ da 
\dz 



I r dx 

R{z) + Q{z) atot - a / —R{x) 
where atot is the NLO total cross section 

/■I dx 

atot^B + aV + a —[R(x)-B]. 
Jo X 



(3.40) 



(3.41) 



As before, we may expand in a to find 

'da\ a 



dz J . z 

naive 



[R{z) + BQ{z)] + Oia") , 



(3.42) 



which shows again that the naive subtraction method suffers from double counting. 

For the modified subtraction method we use eq. (3.38) and eq. (3.39) in eq. (3.20) 
to find 



/ da 
\ dz 



R{z) - BQ{z) 



+ Q{z) [ atot - a dx 



X 



(3.43) 



Expanding in a we now find 



da 

dz 



aR{z) 



(3.44) 



thus, as in the case of the observable y, the MC@NLO prediction at 0{a) coincides 
with the NLO result for all z. The small- ;2 limit of eq. (3.43) is 



lim2; ( — 1 



- *o \dz J 

\ / msub 

From eq. (3.37) we obtain instead 

lim2; I — I 
.-0 \dzj^^ 



aatot ■ 



(3.45) 



aB . 



(3.46) 



Therefore, at z — > the MC@NLO result is equal to that of the standard MC, times 
atot/B, which is the K- factor for the total rate. 
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3.5 Results 

In this section we present toy model results for NLO, MC, and MC@NLO. The 
following parameters have been used in eqs. (3.1), (3.2) and (3.3): 

a = 0.3 , S = 2 , y = 1 . (3.47) 
For the real emission distribution in eq. (3.3) we have taken 

R{x) =B + x{l + x/2 + 20x'^) , (3.48) 
so that the total cross section is 

atot^B + a(v+^^ . (3.49) 

Note that this imphes a large K- factor, Otot/B — 2.34. Finally, for the MC cutoff 
parameter we used 

xq = 0.02 , (3.50) 

and for the maximum energy function in eq. (3.17) we chose the form Xm(x) = 1—x. 
In the left panel of fig. 1 we show the NLO result obtained with eq. (3.13) for the 
observable y defined in cq. (3.22) (solid histogram), compared to the result obtained 
by running the MC (dashed histogram), with Q{x) = 1; the integral of the latter 
result has been normalized to the NLO rate (i.e., each event is given a weight equal to 
(Tfot/N). As a cross check for the numerical implementation against analytical results, 
we also plot as solid hues the functions given in eqs. (3.23) and (3.24) (for the latter, 
B has been substituted with atot- in order to be consistent with the normalization 
of the dashed histogram), which in fact appear to be in perfect agreement with the 
corresponding histograms. The shapes of the NLO and MC results are completely 
different. In fact, the function R{x) has been chosen to introduce a clear difference 
between the NLO and MC results, for a better understanding of the effect of the 
matching between NLO and MC. 

To study the matching of NLO and MC, we consider an MC with a dead zone. 
In our toy model, this corresponds to an MC that cannot emit photons with energy 
X > Xdead (in the numerical simulations we set Xdead = 0.6). We can easily get a toy 
MC with a dead zone by suitably choosing the function Q in eq. (3.14). We first 
introduce a smoothing function of the form 

^(^') = ^2a + c2(l-a;)2/3' ^^-^^^ 

where a, (3 and c are free parameters, and then consider three different implementa- 
tions of the dead zone: 

1. Q{x)^Q{Xdead-x)] (3.52) 

2. Q{x) = Q{xdead - x)G{x/xdead), with Q; = 1, /3 = 1, c = 1; (3.53) 

3. Q{x) = Q{xdead - x)G{x/xdead), with a = 2, /3 = 1, c = 8. (3.54) 
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Figure 1: Left panel: Toy model NLO (solid histogram) and MC results (dashed his- 
togram), overlaid on analytical results of eqs. (3.23) and (3.24), with Q{x) = 1. Right 
panel: MC results, in the case of the dead zone with Xdead = 0.6. Q is chosen as in 
eq. (3.52) (solid), eq. (3.53) (dashed), and eq. (3.54) (dotted); the first bin is normalized 
to unity. 



The results for the pure MC evolution are presented in the right panel of fig. 1; 
the dependence upon Q is actually fairly big, since we made extreme choices in 
eqs. (3.52)-(3.54). It is thus striking that the corresponding results obtained with 
the MC@NLO, shown in the left panel of fig. 2, do not display big differences from 
each other. All three MC@NLO curves are quite close to the pure NLO result (and 
coincide with it for y > Xdead), and differ sizably from it only in the region of small 
X, where the effect of multiple emission is expected to be important.*^ In this region, 
however, the three MC@NLO curves arc close to each other, notwithstanding the 
differences shown in the right panel of fig. 1. We are here in a region where the 
leading logarithmic behaviour of the MC is highly dominant, and thus the impact 
of Q (which is subleading in logarithmic accuracy) is moderate. We see that the 
curve corresponding to eq. (3.52) displays a less regular behaviour than the others 
at y ~ Xdead', however, the effect is very moderate, being at worst of O(a^). 

The results for the inchisive z distribution are presented in the right panel of 
fig. 2. Again, we can see that the three MC@NLO curves are remarkably close to each 
other, in spite of the differences in the implementations of the dead zone. The ratio 
of the cross section in the leftmost bin of the solid histogram over the NLO result for 
the same quantity is 2.23. This value is fairly close to the K-factor (Jtot/B = 2.34, 
that we should find according to eqs. (3.45) and (3.46); the small difference is due to 
fact that this bin is not quite at the hmit z ^ 0. 

^Since y is an exclusive variable, and the total rate is that of the NLO result, the differences 
between the MC@NLO and NLO curves at y > xq are in fact compensated by the region y < xq in 
the total rate. 
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Figure 2: Toy model MC@NLO results; a dead zone is present {xdead = 0.6) and Q is 
as in eq. (3.52) (solid), eq. (3.53) (dashed), or eq. (3.54) (dotted). The results are for the 
exclusive observable y (left panel), and for the fully inclusive observable z (right panel). 
The normalization is absolute, and the dot-dashed histogram is the pure NLO result. 

We can thus conclude that MC@NLO reduces the numerical dependence upon 
subleading logarithms, which is seen when running the MC alone; this is because 
those regions where the subleading terms are becoming numerically important get 
sizable contributions from the hard NLO emissions in the MC@NLO, which are not 
present in the MC, and which eventually become dominant. On the other hand, in 
the soft- dominated regions, the shapes of the MC@NLO curves just reproduce those 
of the MC curves, as we expect. 

4. The QCD case: W+W~ production 

We now face the problem of constructing an MC@NLO in QCD. In order to simplify 
as much as possible all the technical details related to the colour structure of the 
hard process, we consider here the production of W'^W" pairs in hadronic collisions. 
This process also has a fairly simple structure of soft and coUinear singularities. 
However, we argue that the method we propose is actually valid for any kind of hard 
production process. 

We shall use HERWIG [1] as our reference parton shower MC. Other choices are 
of course possible, and our approach is not tailored to work with HERWIG, although 
some of the technical details, reported in the appendices, do depend on the specific 
implementation of the hard process and parton showering in that program. 

In order to define our MC@NLO for W pair production, we rely on the results 
previously presented for the toy model. First we emphasize the main points of the 
toy MC@NLO. Then we make a parallel between the contributions to the toy model 
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cross section, and the contributions to the W+W cross section. Using these results, 
we proceed to construct the QCD MC@NLO, by analogy with the toy MC@NLO. 

4.1 Toy model revisited 

Let us consider again our master equation for the toy model MC@NLO, eq. (3.20). 
The MC@NLO approach consists in running the ordinary MC, with initial conditions 
and weights determined by the NLO calculation. Depending upon the initial condi- 
tions, we distinguish two classes of events, which correspond to Imc{0,Xm{x)) and 
Imc{0, 1). We denote the class corresponding to Imc{0,Xm{x)) by M, which stands 
for "hard", and reminds us that one emission already occurred, due to NLO, prior to 
the MC run. The class corresponding to Imc{0, 1) is denoted by S, which stands for 
"standard" ; this reminds us that this initial condition is the same as for the ordinary 
MC runs. 

It is important to realize that for S events the finite part of the virtual contribu- 
tion, aV, the real counterterm, —aB/x, and the term aBQ{x)/x, serve exclusively 
to fix the normahzation of this class of events, relative to that of the class H. In 
other words, it is irrelevant that the latter two terms explicitly depend upon x: in 
the MC@NLO, the spectrum of the observable O resulting from § events is entirely 
due to MC evolution. In contrast, the NLO hard emissions do contribute to the 
kinematics of H events, and thus help to determine O. 

In order to obtain in QCD a structure analogous to that outlined above, we 
therefore have to collect all the terms of the NLO cross section that do not contribute 
to hard emissions, and use them to compute the normalization of class S events; these 
events result from MC evolutions whose initial conditions are identical to those of an 
ordinary MC. On the other hand, the kinematics of NLO hard emissions determine 
the initial conditions for MC evolutions for IH events. 

Finally, we have already seen in the case of the obscrvables y and z (see also 
app. B.3 for a more formal proof) that at 0{a) the dependence upon Q{x) disappears. 
Roughly speaking, this happens because the contribution due to the Born term B, in 
class S events, cancels at 0{a) with the contribution of the terms ±aBQ{x)/x. We 
observe that, in the toy model, aBQ{x)/x is the 0{a) term in the expansion of the 
MC result (see eqs. (3.24) and (3.37)). We therefore guess that the role of aBQ(x)/x 
in a QCD MC@NLO is played by the 0{as) term in the expansion of the result of 
the ordinary MC, when only the shower is taken into account (i.e., the hadronization 
step is skipped). We shall eventually prove (see app. B.3) that this is in fact the 
case. 

4.2 NLO cross section 

We now deal explicitly with W+W" production in hadronic collisions, and we start 
by writing the NLO cross section. In this paper, we use the results of ref. [22], which 
we rewrite in a form more suited to our present goal; it should be clear that any 
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W"'"W cross section (or at least, one based upon the subtraction method) can be 
cast in the same form. Owing to the factorization theorem, we have 

da{P,,P2) = J2f dx,dx2fi'''\x,)f^'''\x2)daab{xiPi,X2P2), (4.1) 

ab 

for the coUision of hadrons Hi and H2 with momenta Pi and P2 respectively. The 
sum runs over those parton flavours that give partonic subprocesses contributing to 
W pair production at NLO. We symbolically denote these processes as follows: 

qq — > W+W-, (4.2) 
qq — > W+W-g, (4.3) 
qg — > W+W-q, (4.4) 

where q = u, d, . . . and eq. (4.4) includes the corresponding antiquark processes. 
Equation (4.2) contributes to Born and virtual emission terms, whereas eqs. (4.3) 
and (4.4) contribute to real emission terms. 

The information on the hard process in eq. (4.1) is contained in the subtracted 
partonic cross sections, daab (the hat means that counterterms are implicit in the 
expression). We write these cross sections as follows (eq. (3.45) of ref. [22]): 

da^b = da^ + da'--^ + c^aif' + ^^S"^ + da^^ , (4.5) 

where: 

a) da^^^ is the Born contribution. It has 2 — > 2 kinematics, and corresponds to 
eq. (3.1) of the toy model. 

b) da'^^^^ is the contribution of the virtual diagrams, and of all those terms, result- 
ing from various subtractions, that have 2 — > 2 kinematics and are of 0{as). 
This term corresponds to the finite part of eq. (3.2) of the toy model. 

c) da^^i^^ are the contributions of the finite remainders resulting from the subtrac- 
tion of the initial-state collincar divergences (+ and — indicate radiation from 
partons a and b respectively). These terms do not have an analogue in the toy 
model, since in the toy model there are no parton densities. They have 2^3 
kinematics, with the light outgoing parton always coUinear to the beam-line; we 
refer to these configurations as quasi-2 — > 2 kinematics. We stress that da^^^ 
are subtracted quantities, and that the counterterms have 2 — > 2 kinematics. 

d) da^J^^ is the contribution of the real emission diagrams, with its divergences 
suitably subtracted; it has 2 — > 3 kinematics (the counterterms have 2 ^ 2 or 
quasi-2 — > 2 kinematics). The unsubtracted part of this term corresponds to 
eq. (3.3) of the toy model. 
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We point out that all the terms on the r.h.s. of eq. (4.5) are separately finite. Also, 
they do not contain any large logarithms, and so no cancellation of large numbers 
occurs among them. This property is a feature of the subtraction method, and does 
not hold in the slicing method. 

Items a)-d) above achieve the identification between terms contributing to the 
QCD cross section and to the toy model cross section. We now go into more detail, 
by writing the quantities appearing in eq. (4.5) more explicitly . 

We denote by d(f)2{s) and d(l)s{s) the phase spaces of the W"'"W~ pair and of 
the W+W~ pair plus a light parton, respectively; they depend upon the kincmatical 
variables of these particles, and upon s = {pi + P2Y ^ the cm. energy squared of the 
partons that initiate the hard scattering, whose momenta we denote by p\ = XiPi 
and P2 = X2P2- The Born term is 



da'^^{pi,P2) = M'^^{pi,p2)d(f)2{s), 



(4.6) 



where Ai^ab Born matrix element squared, summed over spin and colour degrees 

of freedom (the average factor for initial-state ones is understood), times the fiux 
factor l/(2s). The soft-virtual contribution is 



dcrab\Pi^P2) = 



Qa6A^l?(Pl,P2) +Mir'^(pi,P2) dMs), 



(4.7) 



where Qab is a constant with respect to the kinematics of the partonic process, and 
•^ib "^^^ is the finite part of the virtual contribution. As with any other finite part 
resulting from the cancellation of singularities, its definition is to a certain extent 
arbitrary; however, any redefinition is compensated by suitable changes in the other 
terms contributing to the cross section. The initial-state coUinear terms are: 



where 



dcTab \PUP2) 



2^ 



d(T 



(c±) 
ab 



da 



(c±) 
ab 



ct 



log 



(4.8) 

M^^ {xpi,p2)d(l)2(xs)dx, 

(4.9) 

•^ac (Pl' Xp2)d(t)2{xs)dx. 

(4.10) 



In eq. (4.8), the term da'^^^^^t (the "ct" stands for counterterm) is constructed in 
order to cancel the divergences present in the coUinear "event" contribution da^^^'' |cv; 
these divergences are soft (x ^ 1), as can be seen explicitly in eqs. (4.9) and (4.10). 
The counterterm has a form that is fully specified in a given subtraction scheme (see 
appendix A). In eqs. (4.9) and (4.10), ji is the factorization mass scale (which also 
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enters in the parton densities), and P^^"* is the i*^ term in the e-expansion of the 
Aharelh-Parisi kernel in 4 — 2e dimensions, taken for a; < 1 (i.e., no S{1 
are included). Finally, we write the real emission term as 



ah 



ct 



x) terms 

(4.11) 

(4.12) 



where M.)J is the matrix element squared resulting from the real emission diagrams 

(summed over all spins and colours, and averaged over those relevant to incoming 
partons), times the flux factor. This quantity has soft and initial-state collinear sin- 
gularities, and thus da)^^^ \cv has divergences, which arc cancelled by the counterterms 



da^^^j\ct in eq. (4.11). As in the case of d(J^^^^'\ct• these counterterms are fully speci- 
fied only in a given subtraction scheme. However, one should keep in mind that the 
physical cross section is always strictly independent of the choice of the subtraction 
scheme. 



4.3 Observables at NLO 

As in the case of the toy model, a preliminary step in the construction of the 
MC@NLO is the prediction of an observable at the NLO. We start by defining an 
unintegrated, fully-differential cross section 



dllab{xi,X2) 



ff'\x,)fl"'\x2)daa^{x,P,,X2P2), 



(4.13) 



and analogous quantities, for all the contributions to the partonic cross section that 
appear on the r.h.s. of eq. (4.5). Any observable O can depend, at the NLO, upon 
the kinematical variables of the W+, of the W", and of the fight outgoing parton. We 
denote by 0(3) the observable computed in non-soft and non-collincar kinematics; 
this corresponds to 0{x) in the toy model. We also denote by 0(2) and 0(2) the 
observable computed in 2 — ^ 2 (i.e., soft) and in quasi-2 2 (i.e., collinear non-soft) 
kinematics respectively; these correspond to 0(0) in the toy model. In all cases, 
we understand the variables of the outgoing particles to be given in the laboratory 
frame. We then have 



ab 



dx\dx2d^2, 



0(3)^ 



+ 0(2] 



I, 



+ 0(2)- 



I2 \ d4>2dx 



-{0(2), 0(2)} 



#3 



+ 



ct. 



dE 



ab 



d4>2dx 



dj:?} ^ dT^ 

d(p2 

{c+) 



0(2)- 



1 /rfS 



ab 



T2 \ d4>2dx 



ab 



ct 



d4>2dx 



ct, 



(4.14) 



where in the last term {0(2),0(2)} indicates that the counterterms of the real 
contribution can have both 2^2 and quasi-2 — > 2 kinematics. The factors X2 
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and X2 guarantee that the 2^2 and quasi-2 2 terms are correctly normahzed 
upon integration over the full three-body phase space; the following equations must 
therefore hold 



the integration being performed over the whole phase space. Eq. (4.14) should be 
compared with eq. (3.13). The factors I/X2 and I/T2 ^'^'^ '^'^^ present in the toy 
model, simply because the integral over the phase space dx in that case is equal 
to 1; if this were not the factor analogous to I/I2 would have to be in- 

serted in eq. (3.13). The term c?E^{^/d03|ev —fafbM.^ab corresponds to aR{x)/x in 
eq. (3.13). The term dTi^^^ /d(f)2 —fafbM.^ab corresponds to the Born contribution, 
while the term dr!"^^'' /d(j)2 = fafb{Qab-Mab +•^11''^^^'') corresponds to V, the finite 
part of the virtual contribution. As we said already, the contributions dT,^^^^ have 
no analogues in the toy model, simply because there are no parton densities there. 
Finally, dTiJ^^ /dcpslct represents the counterterms necessary to cancel the divergences 
of the real matrix element; it is thus equivalent to the last term in eq. (3.13), aB/x. 

4.4 Matching NLO and MC: the QCD MC@NLO 

Eq. (4.14) is very similar to eq. (3.13). We could try to follow the same steps that 
brought us from eq. (3.13) to eq. (3.20): replace the observables in eq. (4.14) with 
suitable interface-to-MC's and then modify the equation with a quantity analogous 
to that introduced in eq. (3.21). There is, however, a technical problem, related to 
the initial conditions of class S events. In the toy model, these initial conditions 
are unambiguously determined by letting x — in the corresponding system-plus- 
one-photon configurations. This procedure cannot be automatically generalized to 
QCD, where matrix elements have soft and collinear singularities (corresponding to 
2 and 2 configurations respectively). This may suggest the necessity of considering 
two S classes, defined as the sets of those events resulting from MC evolutions that 
start from 2 and 2 initial conditions respectively. However, the MC cannot handle 2 
initial conditions, since final-state light partons are strictly collinear to the beam-line 
in these cases. 

The solution can be found by observing that, at least as far as the variables of 
the W"*" and are concerned, 2 and 2 only differ by a longitudinal boost. The 
variables of the outgoing light parton arc of no concern here, since they cannot 
contribute to the kinematics of any observables in the configurations 2 and 2, the 
parton being either soft of coUinear to the beam line. The integration over the 
momentum fractions Xi and X2 in eq. (4.14) implicitly contains such a boost. Roughly 
speaking, for fixed Xi and X2 we have 2 7^ 2, but given xi and X2, there exist Xi 




(4.15) 



(4.16) 
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and X2 such that 2{xi) = 2(5j). In practice, we adopt the following procedure.^ 
We treat each term under the integral sign in eq. (4.14) independently (we assume 
some regularization procedure has been adopted in the intermediate steps of the 
computation) , performing a change of integration variables that is in general different 
for each term: 



Zl,Z2), Xi^ x'('\zi,Z2,(l)3) , X2 ^ x''^' (Zi, Z2, (1)3) , (4.17) 



where 03 denotes here the set of independent three-body phase-space variables. The 
functional forms of these changes of variables will be given in app. A. 4. Here we 
are only interested in the fact that eq. (4.17) must be such that 2{xi) — 2{xi). We 
define, for each of the functions T,ab that contribute to eq. (4.14), 



Xn 



d{Zi,Z2) 



-dEi? (4")(zi,Z2,03),4"^(^l,^2,03),03) , (4.18) 



where d{x^^\x^2'^)/d{zi, Z2) is the Jacobian corresponding to eq. (4.17). Then, as we 
shall show in app. A.4, eq. (4.14) can be recast in the following form: 



(O) [ dzidz2d(p3<0{3 

ab ^ 



■(/) 

ab 



+ 0(2) 



T2 



+ 



{sv) 



+ 



1 /rfE 



ic+) 
ab 



d(f)2dx 



+ 



ab 



d(f)2dx 



^2 



ab 



d(t)2dx 



d(t)2 



ab 



d(t)2dx 



ct 



r^(/) 



dH 



ab 



d(f)3 





1 


ct- 





(4.19) 



We stress that the only difTcrcncc between eq. (4.14) and eq. (4.19) is that the latter 
does not contain 0(2) any longer, thanks to the definition of the change of variables, 
which ensures that 2(xj) = 2(0;^). We also stress that no approximation has been 
made in eq. (4.19); it is just a different way of writing the NLO cross section. 

It should be clear that the necessity for the change of variables in eq. (4.17) is 
due to the simultaneous presence in eq. (4.14) of 0(2) and 0(2), which is in turn 
due to the form of the hadronic cross section in eq. (4.1), where the luminosity is 
factorized. One could directly start with an NLO cross section that does not contain 
0(2); clearly, our formulae hold in this case as well, since one just needs to set 

(2^1,^2) = {Xl.^X2). 

It should also be clear that our goal is to obtain an expression for the NLO cross 
section that only contains 3 and 2 configurations. Although we find that we can 
easily achieve this through eq. (4.17), we stress that other methods could be devised. 



^To the best of our knowledge, this method was invented in order to improve the numerical 
stabiUty of the results of ref. [23] . 
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and that our final formula for the MC@NLO will hold, regardless of the procedure 
adopted to eliminate 2 configurations. 

The quantities that contribute to 0(3) and 0(2) in eq. (4.19) correspond to 
those that contribute to 0{x) and 0(0) in the toy model, respectively, according to 

(c±) 

the analogies discussed in sect. 4.2. The only exceptions are the terms da^^ , which 
do not have analogues in the toy model. However, as shown in app. A.4, these terms 
must be treated in the same way as the real matrix element coUinear counterterms. 
This is in fact to be expected: da'fj^^ are strictly related to the coUinear counterterms, 
being finite remainders resulting from the subtraction of the coUinear singularities of 
the real matrix elements. 

Equation (4.19) is therefore strictly analogous to eq. (3.13). This allows us to 
extend to the QCD case the concept of the event classes H and S. In turn, this 
impUes that we need two interface-to-MC's, as in the toy model, which we denote 
by /mc(0, 3) and /mc(0, 2). These are the analogues of /mc(0, Xm{x)) and /mc(0, 1) 
in the toy model. We can thus turn eq. (4.19) into a prescription for a MC@NLO, 
by replacing 0(3) by /mc(0,3), and 0(2) by /mc(0, 2) . In this way, we get the 
analogue of eq. (3.19), that is, a "naive" MC@NLO that suffers from double counting 
and difficulties in the generation of unweighted events. However, precisely as in the 
case of the toy model, we can modify the naive prescription in order to get a QCD 
MC@NLO that implements a modified subtraction, and fulfills the requirements 
given in sect. 2.1. 

To accomplish this task, we need to find the analogue of the quantity aBQ{x)/x 
of the toy model. As discussed in sect. 4.1, what we are looking for is the 0{as) term 
in the expansion of the result obtained by running the ordinary MC; we denote this 
quantity by 



(4.20) 

The specific form of the quantity in eq. (4.20) is dependent upon the implementation; 
different MC's wUl induce different forms of dTiab\yic- On the other hand, (niab\mc 
is independent of the procedure through which we eliminated 2 configurations. The 
form we use in the present paper, and all the technical details concerning its deter- 
mination, are reported in app. A. 5. Here, we merely need to note that dT,ab\mc must 

— (/) 

have the same leading singular behaviour as the matrix element contribution rfE^^ 
when approaching soft or coUinear configurations, in order for the former quantity 
to be used as a local counterterm to the latter. As we shall discuss in app. A. 5, this 
fact is not trivial, and we shall need to modify the 0{as) result obtained from the 
MC. For the time being, we treat d'Lab\Mc in the same way as aBQ{x)/x is treated 
in the toy model (see eq. (3.21)): we subtract and add the quantity 



, , dTiab 



d(f)3 



(4.21) 
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under the integral sign, setting n = 3 in the negative term and n = 2 in the positive 
one. This gives a modified subtraction formula analogous to that in eq. (3.20): 
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(4.22) 



Eq. (4.22) is our proposal for a QCD MC@NLO. In order to use this prescription, we 
need to write it in a form more suited to implementation in a computer code, which 
will be done in the next section. We stress that, as in the case of the toy model, 
the dependence upon O in eq. (4.22) is only formal: the MC@NLO generates events 
without reference to any observable. 

4.5 QCD MC@NLO: practical implementation 

The straightforward implementation of eq. (4.22) results in a computer code which 
is not optimally efficient, in particular in the case of unweighted event generation. 
We therefore implement our MC@NLO through the following procedure. 
We start by computing the following integrals: 
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Note that these two integrals are finite, owing to the properties of dT,ab\ 

(^tot — h + Im: 



MC- 



(4.24) 

Clearly, 
(4.25) 



where atot is the total rate at the NLO level. We also compute the integrals Jh and J§, 
defined as in eqs. (4.23) and (4.24) respectively, except that the integrands are taken 
in absolute value (so that, if the integrands were positive-definite, we would have 
li — Ji). We use a modified version of the adaptive integration routine BASES [24], 
which allows one to compute at the same time the integral of a given function, and 
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the integral of the absolute value of that function. BASES is an evolution of the 
original VEGAS package [25]. 

We then make use of the routine SPRING, which is part of the BASES package. 
Roughly speaking, for any positive-definite function f{x), SPRING returns a set of 
values {xi}2=i distributed according to the function /; the number n must be given 
as an input. Thus, SPRING effectively achieves the generation of n unweighted 
events (with weights equal to 1). In our case, the absolute values of the integrands 
in eqs. (4.23) and (4.24) are given as input functions to SPRING. The numbers of 
events generated in the two cases are (possibly with a roundoff to the next integer) 

Nu = N,,t-r^, Ns = Ntot-r^, (4.26) 

where N^ot is the total number of MC@NLO events that we eventually want to 
generate. It is easy, within SPRING, to keep track of the sign of the integrand 
function, before its absolute value is taken. Thus, SPRING achieves unweighted 
event generation, with the weights equal to +1 or -1. 

We symbolically denote the results of the SPRING runs as follows: 

/h '"^"^ {S.,n,njr}l\, (4.27) 
Is '"^'^ {2.,f.,^f}a; (4.28) 

thus, eqs. (4.27) and (4.28) correspond to events belonging to classes H and § re- 
spectively. For a given i, {3j, fj, wf^''} gives the complete description of a three-body 
event; 3^ represents the kinematics (in the lab frame), fj = {cj, 6j, q} are the flavours 
of the initial- {a-i, hi) and final-state (cj) partons, and w^^'^ ~ ±1 is the weight of the 
event returned by SPRING. The event {2j,fi,w) } is analogous, except that it has 
2 — 2 kinematics (and thus fj = {cj, hi}). 

The events in eqs. (4.27) and (4.28) are the initial conditions for MC evolution. 
In the case of H events, eq. (4.27), the kinematics 3^ is determined by {zi, 4>3). The 
choice of the flavours fj is more involved. One possibility is to flx the pair of flavours 
(a, b) that appears in the sum in eq. (4.23) (and thus to fix f , since c in this case is 
unambiguously determined by a and b), to compute the integral, and to repeat the 
computation for all possible flavour pairs. This procedure is however inconvenient 
from the point of view of computing time. A faster approach is simply to compute 
three integrals of the type in eq. (4.23), for the cases ab = qq, qg and qg summed over 
quark flavours. For each of these integrals, the quark flavours are then determined by 
statistical methods; this can give a wrong assignment of flavours on event-by-event 
basis, but guarantees that the "distribution" in the flavours is returned correctly. 
The reason for not computing a single integral is that subprocesses with qq, qg and 
qg initial states often give contributions of different signs to the cross section, which 
is a source of errors when statistical methods are used. 
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In the case of S events, eq. (4.28), the kinematic configuration 2j is unambigu- 
ously defined as the soft hmit of the corresponding 2—^3 kinematics, determined by 
(2:1, (ps). For the fiavour assignment, we proceed similarly to what was done in the 
case of H events. However, here the relevant flavours are those of the partons enter- 
ing the 2 — > 2 hard process which factorizes in the case of soft or coUinear emissions. 
Thus, we use the Born contribution dT,^^^ to assign the flavour. Notice that this can 
result in errors of 0{as) in flavour assignment. However, the effects induced in this 
way are 0{a^), and thus beyond control. On the other hand, with this prescription 
one achieves at 0{as) the cancellation of the MC terms dHahlMc introduced in the 
modified subtraction, i.e., one avoids double counting. More details on this point are 
given in app. B.3. It should be noticed that, since the Born contribution is always 
positive, there is no need to compute the integral in eq. (4.24) separately for qq, qg 
and qg initial states. 

The above statistical method has the virtue of being easy to generahze to other 
processes. We point out that eqs. (4.25) and (4.26), together with those in the 
remainder of this section, are trivially extended to the case in which more than one 
integral is computed for EI and/or § events. 

The weights in eqs. (4.27) and (4.28) have the following properties: 



and thus, using eq. (4.26), 



-Z^^l -^H, TTl^^i (4.29) 

1 ^ i=i 



Ntot 

We then define the set of events 

{k,, f,, w,}^2r = {3., f., U {2„ f„ }a , (4.31) 

where kj denotes either 3^ or 2j, and 



Wi — w. 



(s,H) 'Jm + Js 



J^^i^'^tot- (4.32) 



We constructed a Fortran code whose output is the set of events {ki,fi,Wi}^2f. 
We point out that the merging of the two event sets on the r.h.s. of eq. (4.31) is 
done randomly; therefore, the probability of finding a S (H) event at an arbitrary 
place in the event list is a constant, equal to N§/Ntot {N^/Ntot)- This set of events 
is passed to HERWIG; there, it replaces the routine that generates the hard process. 
We stress that HERWIG has not been modified; we just added a routine that reads 
the parton configurations coming from the NLO code, and translates them in the 
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standard HERWIG form. Each configuration generates a shower, and all the other 
operations requested by the user. Notice that the unweighting procedure is entirely 
dealt with by the NLO code; HERWIG reads the weights provided by the NLO 
code, and does not discard any events (apart from the few killed in the ordinary 
shower /hadronization process). According to eq. (4.32), in order to get the total 
cross section right, exactly Ntot events must be processed by HERWIG . However, 
should only N^v < Nfot events be processed, the correct rate will be recovered by 
multiplying all the physical results by a factor of Ntot/Nev, since all weights have 
the same absolute value. In other words, any subset of the event set {kj, fj, Wi^^J.'i 
can act as an event set itself, provided that the weights are multiplied by a common 
factor, in order for eq. (4.32) to hold. 

The presence of negative- weight events is a distinctive feature of MC@NLO's. 
Negative weights imply larger statistical errors than in the case of ordinary MC's, 
when the same total number of events is generated. In the HI and S samples of 
eqs. (4.27) and (4.28), the fraction of negative- weight events is 



for i = S, H; clearly, the closer Jj is to /j, the smaller is /j. We can in fact reduce 
the number of negative-weight events in our MC@NLO by varying J§. This quan- 
tity depends upon the form of the count ere vents, or, in a given subtraction scheme, 
upon the free parameters that define the subtraction^. Notice that this is not in con- 
tradiction with the properties of the subtraction method, since J§ is not a physical 
quantity. Thus, we can adjust these free parameters in order to minimize the differ- 
ence Jg — /§; the solution can only be found through numerical methods, because of 
the complexity of the integrand appearing in J§. Further details on this point are 
given in app. A. 3. 

We point out that the integral Jh does not depend upon any free parameters; 
thus, the number of negative-weight IH events is a constant for fixed input (physical) 
parameters in eq. (4.22). We remind the reader that these events will mainly populate 
the hard emission region. 

We remark finally that, owing to the procedure described in sect. 4.4, the 2^2 
kinematic configurations resulting from the integral in eq. (4.24) can all be obtained 
as the soft limit of the 2 — > 3 kinematic configurations identified by {zi, Z2,(f>3)- 
This suggests writing the integration measure in eq. (4.24) with the two-body phase 
space explicitly factorized, d(f)3 — ^02^^03, and integrating over dcp^ before passing the 
integrand function to SPRING. It should be clear that this intermediate integration 
step does not affect our capability of obtaining fully-exclusive final states; however, 
it does have an impact on the number of negative-weight events produced by the 

^These are p and w in the standard subtraction, ( in the C-subtraction: see app. A. 2 and app. A. 3. 




(4.33) 
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MC@NLO. In this paper, we do not exploit this possibihty, leaving the option open 
for future developments. 



4.6 Generalities on matched computations 



We now consider the case of an observable O whose kinematically-allowed range is 
Os < O < Oh- We assume that the region O receives only contributions due 

to hard emissions, and that large logarithms are present in any fixed-order cross 
section for O ~ O5, where a resummation is therefore necessary. In accordance 
with eq. (4.22), we write a differential cross section obtained from the MC@NLO as 
follows: 

da f da\ / da\ f da\ (a oa\ 

Wh~ Waso^ Whs' ^' ' 

The quantity {da / dO)^s is the contribution of all § events, i.e. of the MC evolutions 
-^Mc(0, 2). The quantities {da/dO)ii and {da / dO)i^so are the contributions of the first 
and second terms under the integral sign in eq. (4.22) respectively, and thus both 
are associated with HI events. We expect them to be divergent for O Og, only 
their difference being finite. However, for O — > 0/j, we can assume that {da/dO)^^ 
and {da/dO)^o are separately finite. This is because this region is by hypothesis 
dominated by hard emission, which can only occur through real-emission diagrams; 
this prevents the emitted parton over whose variables we integrate in eq. (4.22) from 
being soft or coUinear. The quantities introduced in eq. (4.34) have the following 
properties: 
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(4.35) 
(4.36) 
(4.37) 
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Eqs. (4.36) and (4.37) imply that, for O — > Oh, we recover the expected result 
da/dO — {da / dO)^^o, which holds up to 0{al) (NNLO) terms. Strictly speaking, 
eq. (4.37) is not necessary for this to happen; {da/dO)^s ^-nd {da / dO)^^^ include the 
leading (and possibly next-to- leading) tower of logarithms log O/O/j, or more pre- 
cisely log[(0 — Os) / {Oh — Os)]-, which become small for O Oh- However, eq. (4.37) 
implies that da/dO ^ {da /dO) ^1^0 also for those values of O which are not too close 
to Oh, but still outside the region where resummation is necessary. The point here is 
the size of the unknown 0{a'^) coefficients. There is in fact no guarantee that these 
coefficients are small; if they are big, the MC@NLO result can differ substantially 
from the NLO one, although the accuracy of the two computations is the formally 
of the same order in perturbation theory. This problem is always encountered when 
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a fixed-order computation is matched with a resummed computation. A possible 
solution (see for example refs. [26, 27]) would be to rewrite eq. (4.34) as follows 
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(4.38) 



where ViO) is a regular function, such that ViOh) — > for O ^ O/j, and ViO) —>■ 1 
for O — > Os, up to terms suppressed by powers of (O — Os)/{Oh — Og). The pertur- 
bative and logarithmic accuracies of eq. (4.38) are the same as those of eq. (4.34); 
however, the function V can be suitably chosen in order to suppress the "MC" con- 
tribution in the region in which the NLO computation can be trusted. Notice that 
the presence of 'D{0) implies that the total rate predicted by eq. (4.38) differs from 
that predicted at the NLO level, even when no cuts are applied. 

The problem with eq. (4.38) is that the function V is specific to the observable 
considered; thus, eq. (4.38) cannot be implemented directly in an MC@NLO. How- 
ever, one could try to devise a solution, similar to that of eq. (4.38), that works 
in the context of MC@NLO's. Clearly, in the hard-emission region the function V 
suppresses the contribution of S events, and that part of H events that corresponds 
to the 0{as) expansion of the MC result. In order to achieve this result without 
any reference to a specific observable, a modification of the showering mechanism 
in the MC is necessary. In particular, one could multiply the Altarelli-Parisi kernel 
appearing in the Sudakov by a function Q{S,^ z) of the showering variables ^ and 
that tends to 1 in the soft and/or coUinear limits, and to zero in the hard emission 
region. It should be noted that this function is completely analogous to the function 
Q{x) that we used in the toy MC, eq. (3.14). Exactly as in the case of the toy model, 
in order to satisfy the requirement that the 0{as) expansion of the MC@NLO result 
be identical to the NLO result, we would need to perform the formal substitution 



in eq. (4.22). Since this procedure would modify the integrand in eqs. (4.23) and (4.24), 
it would also have an impact on the number of negative-weight events generated. In 
this paper, we shall not pursue this issue further, since we prefer to keep the parton 
shower MC unmodified. 

4.7 Results: W+W" and jet observables 

In this section we present results obtained with our W"'"W~ MC@NLO, implemented 
as described in sect. 4.5. We do not seek to give a phenomenological description 
that corresponds to a specific experimental configuration; rather, we wish to show 
the differences between MC@NLO, standard HERWIG MC, and NLO results. For 
this reason, we also present distributions obtained by integrating over the whole 
phase space, i.e., in most cases we do not apply any acceptance cuts. In the case 
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in which cuts are apphed, we only consider the and W~ variables, and not the 
variables of their decay products. We also consider jet distributions, the jets being 
reconstructed by means of the /cx-clustering algorithm [28] . We show results for pp 
collisions at = 14 TeV. We adopt the low-ctg set of MRST99 parton densities [29], 
since this set has a Aqcd value which is rather close to that used as HERWIG default; 
this value of Aqcd (Af^ = 164 MeV) is used for all our MC@NLO, MC, and NLO 
runs. In the case of standard HERWIG MC runs, we give each event (we generate 
unweighted events) a weight equal to o tot I ^ tot-, with Nt^t the total number of events 
generated. For fully exclusive distributions with no cuts applied, this is equivalent to 
normalizing HERWIG results to the total NLO rate Otot- All the MC@NLO and MC 
results (but not, of course, the NLO ones) include the hadronization of the partons 
in the final state. 
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Figure 3: MC@NLO (solid), HERWIG (dashed) and NLO (dotted) results for the rapidity 
of the W"*". HERWIG results have been normalized as explained in the text. 



We present in fig. 3 the rapidity of the W"*", with (right panel) or without (left 
panel) a cut on the transverse momentum p^""" . The sohd, dashed and dotted his- 
tograms show the MC@NLO, MC and NLO results, respectively. This distribution is 
fairly inclusive, and we expect it to be predicted well by NLO QCD in a wide range. 
Prom the figure, we see that the NLO and MC@NLO results are extremely close to 
each other in the whole range considered, whether or not there is a 'p^^ cut. The MC 
result is also similar, but clearly broader than the other two. A similarly inclusive 
variable is presented in fig. 4. The transverse momentum of the W^ is shown with 
(right panel) or without (left panel) a cut on the rapidity y^^ . The previous comment 
applies in this case as well: NLO and MC@NLO basically coincide, whereas the MC 
result is softer than the other two. The same pattern can for example be found for 
the invariant mass of the W pair. We conclude that, for this kind of observable, the 
©(ctg) effects are very small, and NLO and MC@NLO are almost equivalent. This 
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Figure 4: As in fig. 3, for the transverse momentum of the W+. 



also implies that any possible reshuffling of the momenta, due to the hadronization 
phase in MC@NLO, has negligible impact on the colourless W's, as we expect. 

We now turn to the case of more exclusive quantities, such as correlations between 
W+ and W~ variables. In fig. 5 we present p^.^'"^ the modulus of the vector sum 
of the transverse momenta of the two W's. NLO computations cannot predict this 
observable in the region ~ 0, because of a logarithmic divergence for 0; 
on the other hand, NLO is expected to give reliable predictions at large px^^\ The 
MC behaves in the opposite way; thanks to the cascade emission of soft and coUinear 
partons, it can effectively resum the distribution around p^^^ = 0; however, its 
results are not reliable in the large-px*^^ region, which is mainly populated by events 
in which a very hard parton recoils against the W pair. Apart from the unreliability 
of the MC shower approximation to the matrix element in the large-p^,^^^ region, this 
region cannot really be filled efficiently by the shower; this is easily understood by 
observing that, at fixed partonic cm. energy, the large-pJ^*^' region coincides with 
the dead zone (the border of the dead zone in terms of p?^^' can be worked out by 
using the results of app. A. 5). 

This complementary behaviour of NLO and MC approaches can be seen clearly 
in fig. 5, independent of the cuts on the rapidities and transverse momenta of the 
bosons. In the tail, the NLO cross section is much larger than the MC, simply because 
hard emissions are correctly treated only in the former. For p^^^^ — > 0, the difference 
between the two histograms shows the effect of all-order resummation; clearly, no 
meaningful comparison between NLO and data can be attempted in this region. It is 
therefore reassuring that the MC@NLO result interpolates the MC and NLO results 
smoothly. In the small-px^^' region, the shape of the MC@NLO curve is identical 
to that of the MC result. This is evidence of the fact that MC and MC@NLO 
resum large logarithms at the same level of accuracy, as argued in app. B.4. When 
^cyvw) gpQ^g large, the MC@NLO tends to the NLO result, as expected. Again, the 
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hadronization has no significant impact on the W+W system. 

In fig. 6 we present the distribution of the difference between the azimuthal 
scattering angles (i.e., those in the plane transverse to the beam direction) of the 
W+ and W^. This distribution cannot be reliably predicted by fixed-order QCD 
computations in the region A0'^^^ ^ vr; in fact, the NLO prediction diverges loga- 
rithmically for A0*^*' — > IT. This can be seen in the insets of the plots, where the 
variable A^^"^^' has been plotted versus (tt — A0'^^'^^'')/7r on a logarithmic scale, in 
order to visually enhance the region A^'"^^^ ~ tt. We can see that in this region 
the MC@NLO and MC results have identical shapes, as in the case of the observ- 
able p!^^^^ near zero discussed above. The other end of the spectrum, i.e. the tail 
^0(ww) ^ might be considered as analogous to the region of large in that 

this region is dominated, in an NLO computation, by hard single-parton emission. 
However, there are important differences from the point of view of MC treatment. As 
discussed earlier, the region of high p^r^^^ is not filled efficiently by the MC shower. 
On the other hand, the A0'^^^ ~ tail for low pip**' is more easily populated by 
the MC. This is confirmed in the left-hand panel of fig. 6, where we can see that the 
NLO and MC results are quite close to each other at A^^^"^' ~ when no cuts are 
applied. Thus, the region A^'"^^^ ~ receives contributions both from hard parton 
emissions (which dominate the configurations in which the W's have large transverse 
momenta), and from multiple soft or coUincar parton emissions (when the W's have 
small transverse momenta). NLO can treat the former but not the latter, and MC 
can treat the latter but not the former. MC@NLO, on the other hand, is expected to 
handle both emissions correctly, while still avoiding any double counting. When cuts 
on the transverse momenta of the W's are applied, the contribution from multiple 
soft or coUinear emissions becomes less important, as may be seen in the right-hand 
panel of fig. 6. Overall, the MC@NLO results are close to the NLO predictions in the 
/\0(ww) ^ Q shall show later, the differences are compatible with 0{al) 
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Figure 7: MC@NLO (solid) and HERWIG (dashed) results for the transverse energy of 
the hardest jet in each event (left panel), and of fully inclusive jets (right panel). HERWIG 
results have been normalized as explained in the text. 



effects. 

We now turn to the discussion of jet observables. In fig. 7 we present the trans- 
verse energy distribution of the hardest jet of each event (left) , and of the inclusive 
jets (right). In this case, we only compare MC@NLO and MC results; NLO jet re- 
sults are trivial for this process, since there is only one "jet" , which coincides with 
a parton. In the case of the hardest jet, the same discussion as in the case of p^^^^ 
applies; the MC@NLO resums large logarithms at small £'t in the same way as the 
MC does, and can also treat the large-i^^ region, where the MC fails. In the case of 
inclusive jets, the spectrum diverges for Erj, 0, because of the increasing number 
of jets with smaller and smaller E^. In this region, the MC@NLO and MC results 
(rescaled by the K-factor) coincide. We stress that this is actually the same pattern 
that we have already discussed in the toy model: see eqs. (3.45) and (3.46). 
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Figure 8: Scale dependence of MC@NLO (solid band) and NLO (dotted band) results, 
obtained by varying the rcnormalization and factorization scales in the range /xo/2 < 
/Xr,/Xf < See the text for details. 



The NLO results presented so far have been obtained by setting the renormal- 
ization (/Xa) and factorization (/Xp) scales, that we collectively refer to as NLO scales 
(a*nlo), equal to a reference scale //q, defined as follows 



2 1 



(4.40) 



where and are the transverse masses of the W"*" and W respectively. 

The same scale choice has been adopted to compute the integrals given in eqs. (4.23) 
and (4.24), that is, to obtain the sets of H and § events eventually given to the MC 
as initial conditions. 

In general, we refer to the scales used in the computations of eqs. (4.23) and (4.24) 
as MC@NLO scales (//mcbnlo)- We point out that the NLO and MC@NLO scales 
need not coincide. It actually seems more natural to choose different forms for them. 
The MC@NLO scale should in some way match the scale used in the parton shower; 
although any mismatch between the two is formally of higher order than we are 
considering, a choice of the MC@NLO scale motivated by MC considerations should 
result in smaller coefficients for the 0(as) terms in the regions dominated by hard 
emissions. Notice that, when the NLO and MC@NLO scales do not coincide, the 
total rates predicted by NLO and MC@NLO computations are different. 

We shall not pursue the issue of an appropriate MC@NLO scale choice any fur- 
ther in the present paper. We shall limit ourselves to considering the variation of NLO 
and MC@NLO predictions obtained by setting /i^ = = Ho/2 and /Xr = /Xp = 2/io- 
We stress that, while this procedure is a commonly accepted way to obtain a rough 
estimate of the unknown higher-order terms in an NLO computation, it cannot be 
used to guess the MC@NLO scale that best matches the parton shower scale; for 
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this purpose, functional forms other than that shown in eq. (4.40) would need to be 
considered. Still, the procedure can give an idea of the size of the neglected higher- 
order terms in the MC@NLO formalism. The results are presented in fig. 8 for the 
azimuthal correlation A0^^^'; we have chosen this observable since it is the one in 
which the scale effects are largest. The sohd and dotted bands show the variation 
in the MC@NLO and NLO predictions respectively. We observe a sizable reduction 
of the scale dependence when going from NLO to MC@NLO results. Also, in the 
^0(ww) ^ g ^Yie two bands are rather close to each other, and overlap in the 
case that cuts are applied. This is an indication of the fact that NLO and MC@NLO 
results indeed differ by 0(«s) terms, and that the difference is less significant when 
hard cuts are applied, as we should expect. 

An issue related to the choice of scale for the parton showers is the setting of the 
HERWIG parameter EMSCA, which enforces an upper limit on the transverse momenta 
of partons emitted in the showers. The coherence of gluon emission is ensured in 
HERWIG by using an angular variable to perform the shower evolution. In order 
to avoid double-counting of hard processes, it is then necessary to veto emissions 
with transverse momenta greater than the scale of the primary hard process, which 
is specified by EMSCA. If such emissions were allowed, they would constitute the true 
primary hard process. When the veto is applied, the harder emission does not occur 
but the evolution variable is reset to the value selected for the vetoed emission. In 
this way the constraint that the true hard scale is EMSCA is correctly taken into 
account together with QCD coherence in the shower. 

In our case, when running the standard HERWIG MC for W+W~ production, 
we set EMSCA=Y^. In the case of the MC@NLO, there is no reason to have the same 
EMSCA for events belonging to the classes H and §. Wc have ah^cady pointed out that 
the class § must behave in the same way as the ordinary MC, in order to reproduce 
(shape-wise) its results in those regions where a resummation is needed. Thus, for 
S events, we set should also set EMSCA=y^. As far as HI events are concerned, we 
notice that the initial condition for the MC in the toy model can also be interpreted 
as a veto scale; in the case of toy HI events, the initial condition is Xm{x), which is 
imposed to be a monotonic function, with Xm(0) = 1 (see sect. 3.3 and sect. 3.4). A 
possible choice that reproduces this condition in the real case is EMSCA=y^ — 2pT, 
where is the transverse momentum of the final-state light parton, which is always 
present in EI events. Then, in the soft or coUinear limits, this scale tends to that 
adopted in the case of S events. However, such a choice implies that if a very hard 
emission (i.e., Pt — a/s/2) occurs at the NLO, only a very limited possibility is left 
for the MC to emit further, which can lead to problems in the hadronization process. 
Other choices are of course possible. In the toy model, it can be proven that different 
Xm's only affect logarithms beyond the accuracy of the MC. In the real case, we have 
verified that other reasonable choices of EMSCA do not induce visible differences in 
any of the observables studied. 
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Our MC@NLO results are obtained without soft matrix element corrections [6], 
i.e. we set the HERWIG parameter S OFTME=. FALSE. . By doing so, we obtain the 
0{as) HERWIG result reported in app. A. 5, which is crucial in our implementation 
of MC@NLO to avoid double counting. Furthermore, the purpose of soft matrix ele- 
ment corrections is to ensure smooth matching on the boundary between the regions 
dominated by parton showers and by hard emissions, whereas in our formalism there 
is no such hard boundary. In order to have meaningful comparisons with standard 
MC results, we set SOFTME= . FALSE . in that well. We verified that the choice 

of SOFTME has in general only a small impact on the MC results for the observables we 
studied; the only visible effects are in the shapes of distributions in those regions in 
which a resummation is needed: p^^^^ ~ 0, A<p'-^^^ ~ tt, and px^*-* ~ for the hard- 
est jet of the event. A non-negligible effect (~ 10%) is also seen at A^*^^' ~ 0, which 
is partly related to the behaviour at A^^"^^) ~ tt, since the total rate is independent 
of the value of SOFTME. 

We should mention finally that we performed many studies of the stability of 
the MC@NLO results. We studied the dependences of the physical results upon 
the free parameters a and [3 introduced in eqs. (A. 86) and (A. 87), and found them 
completely negligible. This is a rather direct confirmation of the claim in app. A. 5, 
that the angular distribution of soft partons is irrelevant in the determination of 
physical observables. We performed test MC@NLO runs with increased values of the 
showering cutoffs VQCUT and VGCUT, and found no significant variations with respect 
to the results presented here; this supports the claim that MC and MC@NLO have 
the same dependences upon MC showering cutoffs. Also, we find that the parton 
densities used for MC backward evolution do not strongly influence the MC@NLO 
results for this production process. The results obtained using low-ctg MRST99 in 
the MC evolution are statistically compatible with those obtained using MRST-LO 
(low-CKs MRST99 being used in both cases for the NLO part). 

5. Other processes 

The extension of MC@NLO to other processes, at least those with unambiguous 
colour flow, seems straightforward in principle. The key equation, (4.22), is formally 
valid for any kind of singularities. For flnal-state coUinear singularities, the choice of 
the Bjorken x's and the flavour assignment are trivial: one just has to give the MC 
the lower-order process that preceded the collinear splitting. 

In processes where there are several possible colour flows, the procedure adopted 
in the HERWIG MC is first to select the kinematics according to the full differential 
cross section, including interference between colour flows, and then to assign a colour 
flow according to the relative contributions to the cross section neglecting interfer- 
ence (i.e. in the large- A^c limit) [30]. This procedure can easily be extended to the 
MC@NLO case. We point out that the colour flows can be determined independently 
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for S and HI events. The correct NLO results will still be recovered upon expansion in 
tts, since the cancellation of the MC terms (4.21) inserted in the modified subtraction 
formula, eq. (4.22), occurs on statistical basis, and not event by event. In order to 
construct the MC-like cross sections (4.20), the technique presented in this paper can 
be adopted: namely, the pure MC result is smoothly matched onto that predicted 
by perturbative QCD in the soft hmit. In fact, independently of the choice of the 
showering variables, the former cross section is in general locally different from the 
latter, since it is only in the former that the radiation function is integrated over an 
azimuthal angle, in order to implement colour coherence in a MC-friendly manner. 
For processes involving emission from gluons, there arc also azimuthal correlations 
in the collinear limit, which integrate to zero and are therefore not included in MC's 
but cannot be neglected locally. 

For some processes, such as QCD jet production, the total cross section atot is 
not defined without a cut; for example one may require the total amount of trans- 
verse energy to be bigger than a given quantity. The formulation of the MC@NLO 
described in this paper should remain valid in the presence of such a cut. 

In processes involving incoming photons, pointlike and hadronic components 
needs to be defined in a given subtraction scheme. Within the NLO approximation, a 
change of scheme gives rise to terms that are of higher order. However, for consistency 
the same scheme should be adopted for the generation of the NLO cross section and 
the backward evolution of the parton shower in the MC. We also point out that, in 
order for the procedure described in this paper to be applicable, the photons should 
have a continuous energy spectrum, such as that from an incoming electron beam, 
in which the photon energies are distributed according to the Weizsacker-WiUiams 
function. 

6. Conclusions and outlook 

We have presented a practical proposal for constructing an event generator which, we 
believe, displays the main advantages of the NLO and MC approaches. For brevity 
we call this MC@NLO. A key ingredient is the use of a modified subtraction method 
for deahng with infrared and collinear singularities. The method has been illustrated 
using a toy model and implemented in full for the process of hadroproduction of W 
boson pairs. The results there look encouraging (although of course comparisons 
with data are for the future): all the obscrvables investigated show the expected 
behaviour in regions where the NLO or MC results should be most reliable, with a 
smooth transition between hard and soft / collinear regions. 

The method generates a fraction of events with negative weights, but in the 
cases we have studied the fraction is small and unweighted (signed) event genera- 
tion remains rather efficient. One can reduce the number of negative weights for a 
given process by tuning the parameters specific to that process (these parameters 
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have nothing to do with the global MC parameters, which are tuned once and for 
all by comparison to data). We shall consider in future upgrades of the implementa- 
tion of our method the strategies outlined in sect. 4.5 and in sect. 4.6, which should 
reduce the number of S and EI negative-weight events respectively, using an inter- 
mediate integration step, and a modification of the showering algorithm. One could 
also consider more exotic possibilities, such as defining "superevents" by combining 
configurations that are close together in three-body phase space. 

The implementation presented here is based on the HERWIG MC generator, but 
there appears no reason why other parton shower generators such as PYTHIA should 
not be used. The main points to be reconsidered when changing the MC genera- 
tor will be the dead zones for hard emission, the matching of local counterterms 
in the soft region, and the relation between the 2- and 3-body configurations con- 
nected by the backward evolution in the showering initiated by initial-state partons 
(corresponding to the event projection that we used in this paper). 
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A. W+W~ cross sections 

In this appendix, we give all the details relevant to cross sections at 0(0;°) and 0{as). 
Although we study explicitly the case of W+W" production, some of the results are 
based on fairly general ideas, which can be applied to any other hard production 
process. 

A.l Kinematics 

We start by defining some kinematic variables for W+W~ production. The momenta 
in the 2 — > 3 partonic processes are assigned as follows: 

a{pi) + b{p2) c{k) + W+W-(A;ww) ■ (A.l) 
Following ref . [22] , we introduce the quantities 

x = M^^/s, y = cose, (A.2) 

where s — (pi +^2)^, -^ww is the invariant mass of the W pair, and 9 is the scattering 
angle of the outgoing light parton in the cm. frame of the colliding partons. In this 
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frame, 1 — a; is the energy fraction of the outgoing hght parton (in units of 
These quantities unambiguously identify the soft and initial-state coUinear regions 
of the three-body phase space. In particular 

X = 1 =^ soft (k° = 0) , (A.3) 
y — 1 =^ coUinear {k \\ pi) , (A.4) 

— * 

y — —1 ^ coUinear {k \\ P2) . (^-5) 
The ranges allowed for x and y are 

p<x <1, -l<y <1, (A.6) 

where p = im^/s. Defining 

= p,^^/^^, (A.7) 

the two-body and three-body phase spaces are 

d02(s) = -^dcos6'w , (A.8) 

iDTT 

d(f)2{xs) — -^^dcos9^/v : (^-9) 

iDTT 

d^Psis) = — —^{1 — x)dx dy d cos 6yv dip , (A. 10) 

where and (p are the polar and azimuthal angles of the W+ in the W+W~ cm. 
frame. The normalization in eq. (A. 10) is such that both 9^ and (f range between 
and TT. 

A. 2 Standard subtraction 

Here we give the subtraction terms using the definition of ref. [22]. The quantity Qab 
appearing in eq. (4.7) reads 

Q.,= ^(log^(6 + 161og^) + 321og^;g-^7r2) ' ^^'^^^ 
and Qab — for {a, b} ^ {g', q\. We have introduced a free parameter p such that 

p<p<l, /3=0^. (A.12) 



The coUinear terms of eq. (4.8) read 



, UJS ( \ \ _ /logfl - x\ 
log 7T-^ + 2 



1\xr \\ — X ) - \ 1 — X 



-P^\x)'^M'S{xpi,P2)dcl>2Mdx, (A.13) 
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and analogously for da^^ . The real emission contribution of eq. (4.11) reads 

x{l-x){l-y^)M^'\p,,P2)dMs) . (A.14) 

The distributions introduced in eqs. (A. 13) and (A.14) define the subtraction terms 
as follows: 

p.,(.,(_i_)^ 

where < < 2, and h and g arc generic regular functions. It is clear that a variation 
of the (arbitrary) parameters p and ui amounts to a finite redefinition of the soft and 
coUinear counterterms respectively. The cross section is strictly independent of the 
values of p and cu in their whole range (i.e., these parameters do not need to be small, 
as opposed to the slicing method) . 

A. 3 ^"Subtraction 

We give here an alternative form of the subtraction terms, which has been used in 
our implementation of the MC@NLO. Although the results of the MC@NLO do not 
depend on the particular form adopted for the subtraction terms, we find that the 
form introduced in this section gives a superior numerical stability, in a sense that, in 
the context of our MC@NLO, allows to reduce the number of negative- weight events. 
See sect. 4.5 for more details. We define 

V{x,y) = {l-xf{l-y^), (A.18) 

where x and y are defined in eq. (A. 2). The function V has a simple physical meaning; 
in fact 

p^^'^'^^Vn^), (A.19) 

where p!^^^ is the transverse momentum of the W pair, which coincides with that 
of the outgoing light parton at NLO. The ^^subtraction basically amounts to sub- 
tracting the counterterms when the condition 

r{x,y) < C (A.20) 



^ , h(x) - e(x - p)h(i) 

ax -^—^ ^ — — 

1-x 



(A.15) 



^ , h(x) - Q(x - pWl) ^ , 

dx^^ ^ ^' ^ ' log(l -x), (A.16) 

\ — X 
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is met, ( being a numerical parameter, < C < 1. More precisely, we define 



' ^^ hix)-e{C-V{x,y))h{l) 
^ 1-x 



(A.21) 



dx h{x) 



dygiy) 



F{x) 
1-x. 
1 



rfx— ^^^^ \ JJJ y > F{x), (A.22) 



-1 



1-x 

g{y)-m-n^,y)mTy)9{Ti) 

l±y 



, (A.23) 



where y in eqs. (A.21), (A.22), and a; in eq. (A.23), are held fixed. We have 



qq 



in 



log— (6 + 161og/3)+32V/?--7r^ + :^,e(^) , 



2tt 



(A.24) 



log 



2/x^ VI 



+ 



r 



logjl - J^,{x)) 
1 — X 



+ 2 



log(l — x] 
1-x 



V J 



1 - x)P^^\x) - pW(x) \M^^{xp,,p2)d<t>2{xs)dx,{A.2b) 



^(_1_)[(_1 



+ 



l + y 



X (l-a;)(l-|/2)A<(^)(Pi,P2)#3(s) 



where 



yc<i-p, 

v^>l-p, 



(A.26) 



(A.27) 



and 



J^s = 4 log -I log - 2 log log(4(l - f)) 
P 1-1/ 1 - ?/ 



+4Li2 



1+y 



4Li, 



1-^ 



Tc{x) 



1-c/ii-xy x<i-vc, 

a; > 1 - VC- 



(A.28) 
(A.29) 



We point out that the subtraction scheme defined in this appendix is rather general. 
In fact, the condition that determines the subtraction, eq. (A. 20), applies to the 
variables of the parton that recoils against the system, regardless of the nature of 
the system itself. 

We can now go back to the point raised in sect. 4.5, concerning the number of 
negative-weight events in our MC@NLO. As mention there, this number is a (very 
complicated) function of the free parameters used to define the subtractions terms: p 
and cu for standard subtraction, C for ^-subtr action. However, for a crude estimate of 
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their best values, we can assume that the partonic cross sections do not depend upon 
them, and thus that the parameter dependence of J§ is only due to the size of the 
phase-space regions in which the subtraction of the counterterms is performed (this 
amounts to saying that the only parameter dependence taken into account is that 
due to the © functions in eqs. (A.15)-(A.17) for the standard subtraction, and in 
eqs. (A.21)-(A.23) for the ^-subtraction) . In general, there is a cancellation between 
the first two terms in eq. (4.24); to a certain extent, the term — dE^{^/d03|ct can 
be tuned in order to cancel with dTiab/d(t>z\Mc- This explains why C-subtraction is 
a better choice than standard subtraction in the present case: the region in which 
the counterterms are subtracted, eq. (A. 20), is a better match to the dead zone, 
eq. (A. 69), compared with the corresponding region for the standard subtraction. 

In this paper we do not try to find an optimal value for the parameter C,. The 
crude estimate mentioned above gives C, = 0.275; the fraction of negative-weight 
events that we get by using this value is 12%. This fraction appears to be remarkably 
stable against variations of C in the region 0.05 < C i$ 

A. 4 Choice of Bjorken a;'s 

Here we give the form of the change of variables in eq. (4.17). We stress that this is a 
convenient procedure for getting rid of the 2 configurations that appear in eq. (4.14), 
which also has the virtue of being useful for defining the MC-hke terms of eq. (4.21), 
as we shall show in app. A. 5. These two problems could in principle be solved 
independently, and by different means, as it may be necessary in other MC imple- 
mentations; however, our master formula, eq. (4.22), would still hold. We denote 



the three-body configuration that is obtained by using the three-body phase-space 
variables 03, and the momentum fractions Xi and X2; the momenta in eq. (A. 30) are 
understood in the hadronic cm. frame. Furthermore, we denote by 



the two-body configuration that is obtained by setting to zero the energy of the 

(s) 

final-state light parton; means "soft limit" of ^3. In the parametrization given 
in eq. (A. 10), this limit can be obtained by letting x — >■ 1 (see eq. (A. 3)). It has 
to be pointed out that any 2 configuration appearing in eq. (4.14) can be viewed as 
the soft limit of a 3 configuration; thus, eq. (A. 31) is fully general. Analogously, we 
denote by 



by 



3{xi,X2,4>3) 



(A.30) 



2,(xi, X2, 03 ) = S{Xi,X2, 03 ) 



(A.31) 



2c+{xi,X2,(t)3^^) = 3{xi,X2,(t)3^^) 
2e_(xi,X2,0^''~'') = 3(xi,X2,(j)'^~^) 



(A.32) 
(A.33) 
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the three-body configurations that are obtained when the final-state hght parton 
becomes colhnear to the incoming parton coming from the left (pi \\ k) and from the 
right {p2 II ^) respectively; these limits correspond to y ^ 1 and y — > — 1 respectively 
(see eqs. (A. 4) and (A. 5)). 

We now look for changes of variables of the type given in eq. (4.17). There, 
the possibility exists to have a different form for each of the terms contributing to 
eq. (4.14). However, since we are looking for a solution to the equation 2(xj) = 
2(xi), we shall actually have to find a change of variables (4.17) for each kinematical 
configuration 2 and 2. Since only kinematics matters here, it follows that the collinear 
counterterms to dT!'f^\cy. and the terms rfS^^"^^ |cv, must be treated exactly in the same 
way in eq. (4.14), since they all have 2 kinematics. Similarly, the Born term dTj~^l, 
the soft-virtual term rfS^j , the soft counterterm to |cv, and the counterterms 
(iS^^ let, are treated in the same way, since they have 2 kinematics. In summary, 
our task amounts to finding the three changes of variables 

xi ^ x^f^ {zi, Z2Az) , 2:2 = ^2,03) , q; = s,c-F,c-, (A.34) 

such that, at fixed Zi, z^ and ^3, 

1s{x^t\zx, Z2, 03), 4'^(^l> ^2, 03), (t>Z^) 

= 2,+ (a;^Hzi, Z2, <t>,),x^t\zu z^, <f>s), 0?+^) 

= 2c_(xf~'*(^l,^2,03),4''~''(^l,^2,03),0f~'') , (A.35) 

where the equality holds for the and variables only. It is then apparent 
that eqs. (4.14) and (4.19) are identical; in fact, if 2 = 2 for the W"*" and W~ 
variables, then 0{2) — 0(2), since the outgoing light parton cannot contribute to 
the kinematics of any infrared-safe observables, being either soft or collinear to the 
beam line. 

The requirement 2 = 2 amounts to two conditions, since in both configurations 
the transverse momentum of the W pair is equal to zero, and only two components of 
the four-momentum of the pair are left to be constrained. Notice that it is sufficient 
to constrain the pair momentum, since the internal variables of the pair are held fixed 
in eq. (A.35). Instead of imposing eq. (A.35) directly, it turns out to be convenient 
to follow the strategy developed in ref. [23] (where the procedure described here is 
called "event projection"), where 2 and 2 are related to 3. This is done by imposing 

Oi(3(^i,;^2,03)) = 0,{2,{x['\xi'\4'^)), (A.36) 
02(3(^1,^2,03)) = 02{2,{x['\xi'\4'^)), (A.37) 

and analogous equations for c+ and c— ; Oi and O2 are two arbitrary observables, 
defined with the variables of the W+W~ pair; therefore, the conditions in eq. (A.35) 
are effectively implemented by relating 2 and 2 to 3. It must be stressed that 
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eqs. (A. 36) and (A. 37) are imposed in the hadronic cm. frame for those observables 
that are not invariant under fongitudinal boosts. 

Notice that this procedure can be apphed to any production process; in the 
general case, the observables Oi and O2 are defined with the variables of the system 
that recoils against the hard parton emitted by real diagrams at the NLO level. In 
the present case, one obvious choice is 

Oi = Ml^ , O2 = yww , (A.38) 

where and |/ww are the invariant mass squared and the rapidity of the pair 

respectively. These two observables are easily expressed in terms of the momen- 
tum fractions xi,2, and of the phase space variables in the parametrization given in 
eq. (A.IO): 

Ml^{3{zi,Z2,(l)3)) = Ml^{zi,Z2,x)^xziZ2S, (A.39) 

1 Zi 

yww(3(2;i, Z2, 03)) = Vwwizi, Z2, X, y) = log^ + - log — , (A.40) 

2 Z2 

where S is the hadronic cm. energy squared, and 



]j2-(l-x){l-y) ^ ' 

By imposing the conditions in eqs. (A. 36) and (A. 37), and using the explicit forms 
for and yww given in eqs. (A.39) and (A.40), we get the following: 

x?^z,^V^, xi'^^Z2^, (A.42) 

x^^^^ =x['^/x, xf+)=4^), (A.43) 

xt^=x['\ xt^ = x^^/x. {AAA) 

It should be stressed that x^f^ and x^^^ as given in eqs. (A.42), (A.43), and (A. 44) 
do not need to be less than 1 any longer; thus, the parton densities contained in 
have to be defined as equal to zero if their arguments are larger than 1. 

Although the choice made in eq. (A.38) is arbitrary, we shall show later that, 
in the context of our MC@NLO, it is actually dictated by the procedure used by 
the MC to generate the Bjorken x's. In this procedure, the MC typically constrains 
either the invariant mass and the rapidity of the pair, or the invariant mass and the 
longitudinal momentum of the pair. In the current version of HERWIG the latter 
choice is adopted. We therefore also consider 

Oi = Ml^ , O2 = /cww , (A.45) 
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where fcww is the longitudinal momentum of the pair. With the parametrization of 
eq. (A. 10), we get 

k^w{^{zi,Z2,(j)3)) = k^w{zi,Z2,x,y) = — — [(1 - x)y{zi + Z2) + {I + x){z2 - Zi)] . 

(A.46) 

By imposing again eqs. (A. 36) and (A. 37), we now obtain 

x'f^ - ^ (VB^ + AA - , x^'^ = ^ ( + AA + , (A.47) 

a;S'^+) =a;S^)/a;, +) = , (A.48) 

x^r^ = x^\ x^r^ = x^'^/x, (A.49) 

where 

A^xziZ2, (A.50) 

^ = ^ [(1 - ^)y{^i + ^2) + (1 + x)(-Z2 - zi)] . (A.51) 

Note that eq. (A.48) is identical to eq. (A. 43), and eq. (A.49) is identical to eq. (A. 44). 
This gives a consistent picture when the results given above are used in the MC@NLO, 
eq. (4.22). In fact, event projection has been introduced in the context of an 
NLO computation, and the constraints have been imposed on final-state quanti- 
ties only. However, in an MC@NLO the initial-state configuration of the hard 
process is relevant as well, since it defines the initial conditions for the MC evo- 
lutions implicit in 7mc(0, 3) and /mc(0,2). Eqs. (A.42), (A.43), and (A.44) [and 
eqs. (A.47), (A.48), and (A.49)] all enter the definition of § events, and since in 
general xf^ 7^ x^^^'^ 7^ a;-^ \ there is an ambiguity for the prescription of the mo- 
menta of the initial state partons that initiate the backward evolution. However, as 
explained in sect. 4.5, the initial conditions for class S events are the 2 — > 2 hard 
processes; the initial states of these processes coincide, by definition, with those of 
the 2^2 contributions to § events, but do not coincide in general with those of the 
quasi-2 2 contributions. If we denote the momentum fractions of the initial-state 
partons entering the 2^2 processes as x'-^ and x^, owing to the fact that 1 — x is 
the energy fraction of the outgoing light parton (see app. A.l), we have 

^(«)' _ ^(s) ^(«)' _ As) /A Kn\ 

JU — Jb J J-"2 — '2 "J yl\..xj ^ J 

(c+)' ^ _ (s) (c+y ^ (c+) _ (s) / . 

2 — Jb Jb — ju-^ ^ — 2 — 2 ' V"^* / 

(c-)' ^ (c-) _ is) (c-y ^ (c-) _ (s) / . 

Thus, xf"* = xf~^^ = \ and there is no ambiguity in the momenta of the initial 
state partons. It should be noted that an ambiguity here would in any case affect our 
results at 0{a'^), since it would only be due to quasi-2 — > 2 contributions; however. 
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the absence of any ambiguity fits nicely in a picture where MC and NLO give different 
descriptions of the same physics. 

In appendix A. 5, we shall need to consider the soft and coUinear limits of x\^^ 
and Xj-'^^^ By direct computation, we obtain the following 

X2 — >■ ^2 , 

x^l =^ <( ^ , 4'^) ^ Z2 , (A.55) 

(c-) 
X2 Z2 , 

X2 — >■ ^2 , 

(c+) 
X2 Z2 , 

X2 — ^ XZ2 , 
(c-) 

X2 Z2 , 




(A.56) 
(A.57) 



(c— ) fc+) 

where we did not consider the y ^ 1 hmit of xi and the — > — 1 limit of a;^ , 
since these limits are not relevant in the computation of singular cross sections. We 
stress that eqs. (A.55)-(A.57) hold regardless of whether eq. (A. 38) or eq. (A. 45) has 
been used. This is a manifestation of the factorization properties of QCD, as will 
become clear in app. A. 5. 

A. 5 Dead zone and subtraction term 

Here we define the quantity dT^ablMc, introduced in eq. (4.20). We make use of some 
of the results of ref. [8]; although ref. [8] concerns only single vector boson production, 

the kinematics carry over to the present case if we simply replace the momentum of 
the vector boson by that of the pair. According to the discussion of sect. 4.1, we have 
to consider the cross section that we would obtain by keeping the 0{as) term in the 
expansion of the result of the MC, which we denote by da\j,ic- Roughly speaking, 
da\Mc is equal to the Born cross section, times the 0{as) term in the expansion of 
the Sudakov form factor. More precisely 



da 



= V dxi'^'>dxt''^^daf^{x^r'>p^, x^r^P2) 

MC ^ — ' ZTT 

abc 



X 



f±^e(4 - cM^\z^)e^\xr'/z^)fi'''Vr') 

?+ z+ 



dE- dz. 
+ - 



Z- 

Equation (A. 58) has been obtained directly from the QCD Sudakov form factor; in 
doing so, we have eliminated any cutoffs, precisely as in the case of the toy model. 
This impfies that eq. (A. 58) has soft and coUinear divergences, which is what we want, 
since from it we derive dT^ablMc, namely a local counterterm to the real emission 
matrix element. The two terms on the r.h.s. account for the emissions from the 
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incoming partons. If strongly-interacting partons were present in the final state at 
the Born level, corresponding terms should be added. The variables ^± and z± are 
those used by the HERWIG MC [4] in the parton showering process. The functions 
define the region where an emission can occur during the shower; therefore, from 
eq. (A. 58) we can easily read the form of the dead zone, which is by definition the 
region where no emission can occur. We have 

Dead zone : > , > , (A. 59) 

where both conditions must be fulfilled. Finally, the l/z± factor that appears in the 
arguments of the parton densities in eq. (A. 58) is due to the backward evolution in 
the shower from initial-state partons. In the case of emission from final-state partons, 
the arguments of the parton densities would simply read x^^'^'^ and x^2^'^\ 

The values of x''^'^^ and x'^^'^^ are determined by the MC (through energy- 
momentum conservation) at the Born level, after the generation of the 2 configu- 
ration which is used as initial condition for the shower. At the end of the shower, 
the momenta of the final-state particles are boosted, the boost being determined by 
imposing a relation between the variables of the W pair (in general, of the system 
produced at the Born level) before and after the shower. It follows that, after the 
shower, x^^^^"* and x^^'^'^ have still the values determined by the 2 configuration; how- 
ever, after the shower x^^" ' and x^2^'^'^ cannot be interpreted any longer as Bjorken 
x's associated with the incoming partons. 

This procedure is identical (when we consider the 0{as) result, that is, the 
emission of a single parton) to the one we defined for event projection in app. A. 4. 
In other words, the 3 configuration imphcit in the l.h.s. of eq. (A. 58) is related to 
the 2 configuration that initiates the shower by means of constraints such as those 
in eqs. (A. 36) and (A. 37).^ It follows that, if we label the 3 configurations according 
to the conventions of this paper, 3 = 3(^1, 2^2, 03 ), we have 

xS"^) = , 4"°) = . (A.60) 

The explicit form of x^^'^ depends on the observables Oj used by the MC to imple- 
ment eqs. (A. 36) and (A. 37). A priori, these need not coincide with those given in 
eq. (A. 38) or eq. (A. 45). In practice, HERWIG uses the variables in eq. (A. 45), and 
therefore the results in eq. (A. 47) give the explicit form of x^^'^'' and x^^'^\ The 
relations between 2 and zi, Z2 allow us to recast eq. (A. 58) in the following 

form 



da 



= dzidz2dEab , (A.61) 

^ ' MC 

ab 



''This is the reason for introducing event projection in app. A. 4, instead of just imposing 
eq. (A. 35). In a completely general case, one should map the three-body configurations obtained 
from the MC onto the 3 configurations of the NLO cross section, thus defining a dTiab\Mc that acts 
as a local counterterm in eq. (4.22). 



48 



where 



dE 



ab 



MC Z+ 0[Zi,Z2) 



MC 



+ f ^l^/f ^n-S^Vr^l-rV-)^-;. (A.62) 

Z- (j[Zi,Z2) MC 



is the quantity that the appears in eq. (4.22), and 



= £ ^^^+0(4 - ^^)P^\^^)daT{x['^Pu x^^P,) , (A.63) 
= ^ |^cZ._e(.! - A, 4^^P2) , (A.64) 



where the sum over the parton label c is now understood. 

In order to prove that dYlab\mc in eq. (A.62) can act as a local counterterm in 
eq. (4.22), we have to express it in terms of the variables we use to describe the NLO 
cross section. This can be easily done as explained in ref. [8] : we relate the variables 
^± and z± used in the HERWIG parton showering to Lorcntz invariant quantities, 
which we then express in terms of the phase-space variables of eq. (A. 10). Following 
ref. [22], we define 

tk = (Pi - kf = -|(1 - x){l - y) , (A.65) 
Uk = (P2 - kf = -^(1 - x){l + y) , (A.66) 

where the momenta are assigned as in eq. (A.l), and x and y have been defined in 
eq. (A.2). We have [8] 

tk^-^—4^i+Ml^, (A.67) 



tkUk (1 - z+) 



z: 
2zl 



\2 



e+(2-e+)M^^, (A.68) 



for a shower initiated by parton a. A shower initiated by the other incoming parton 
(fe) is described in terms of and z_, which are in turn related to tj. and by 
the same functional forms given in eqs. (A.67) and (A.68), with the formal exchange 
ik ^ Uk- We can now substitute eqs. (A.2), (A.65), and (A.66) into eqs. (A.67) 
and (A.68), to find C±{x, y) and z±{x, y), whose explicit form we do not report here, 
since it is not needed in the following. 

We can now find the dead zone boundary in the {x, y) plane; from eq. (A. 59) we 
obtain 

Dead zone : p < x < x^z , \y\ < Yoz{x) , (A. 69) 
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where both conditions are met simultaneously, and 

= (A.70) 

yoz(x) = 1 - (3 - VT+s^) . (A.71) 

We now consider eqs. (A. 63) and (A. 64). Using eq. (A. 10), we obtain 



dcj>: 



'3 



(A.72) 



d4>3 



MC s (l-x)^- d{x,y) 

(A. 73) 



where pi = ZiPi (according to the definition given in eq. (4.13)); therefore, s = Z1Z2S. 

We can now consider the soft and collincar limits of dHablMc given in eq. (A. 62). 
We observe that, owing to the procedure adopted by the MC to determine x^^'^'' and 
x^'^\ the luminosities that appear in eq. (A. 62) have non-trivial soft and coUinear 
limits. We make use of the fact, obtained by direct computation, that 

lim z± — 1 , lim z+ — x , lim z^ — x . (A. 74) 

Therefore, using eqs. (A.55)-(A.57), we obtain 

= hm/f ^)(4^V^+)/f^^(4^^) = lim /r^H4^Vf^^(4^V^-) 

= /f ^n^i)/rH^2) , (A.75) 

where the coUinear limits have been considered only for the luminosities attached 
to short- distance cross sections which are singular in those limits. In the soft and 
coUinear limits, dHabluc wiU have to cancel the singularities of dulf,^ and dH^^^K 
Owing to the event projection procedure described in app. A. 4, the luminosities 
attached to the latter two terms behave as follows in the soft and coUinear limits: 

a;— >1 >1 

= lun f("^\xt^)ff '^\xt^) = fr^\z,)ff^\z,) , (A.76) 

where we used eqs. (A.55)-(A.57). Eqs. (A.75) exactly match eqs. (A.76). It is also 
straightforward to prove that 

Z± d[Zi,Z2) 
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in the soft and collinear limits, and that the corresponding Jacobians in the NLO 
cross section have the same behaviour. This imphes that the properties of dJ^ahlmc 
as a local counterterm are determined by the behaviour of its short-distance parts, 
given in eqs. (A. 72) and (A. 73). 

When considering eqs. (A. 72) and (A. 73), we notice that these formulae are 
valid for any parton flavours a and b. However, in our case the soft and collinear 
limits are all non-trivial only in the case a — q, b — q, with c — g;we shall therefore 
explicitly discuss only this case. We expect the forms of the leading soft and collinear 
singularities of the cross sections given in eqs. (A. 72) and (A. 73) to be dictated by 
QCD factorization theorem. In other words, we expect the following equations to 
hold 



dcTqg{pi,P2) 



+ 



87rQ;<jCE 



647r«<, 



Pi ■P2 



Pi-k p2-k 



-Mg(Pl,P2) 



-MS?(Pi,P2),(A.78) 



s (l-a;)2(l-|/2) 



d03 



dagg{Pl,P2) 



Pi ■ k 



P^^{x)M^S{xp,,P2) 



dh 



s (1 

-1 47ra 



X) 



y) 



Mfhxpi,p2) , 



(A.79) 



MC 



P2-k^'"' 



:M''q^{puxp2). 



. (l-x)(l + y)'"^^-^'-^^-- ^^-^^^ 

Note that in the case of the soft hmit we sum the short distance cross sections, thanks 
to eq. (A. 75). Explicit computations of the collinear limits y — > ±1 in eqs. (A. 72) 
and (A. 73) give eqs. (A.79) and (A. 80). On the other hand, in the soft limit x ^ 1 
we find the following: 



d(Tqq{Pl,P2) 



d(f)3 



+ 



d(Tqq{Pl,P2) 



MC 



d(f)3 



x^i 64:TTas^ . (b) 



C,M':^{p,,P2) 



MC 



X 



2e(y + l/3) 



+ 



2e(-y + l/3) 



,(l-x)2(l-y)(3 + y) ' (l-x)2(l + y)(3-y); " ^^'^^^ 

It appears therefore that the choice of the showering variables in HERWIG violates 
the factorization properties of QCD cross sections. However, the situation is not as 
bad as it may seem. In fact, the kernel in eq. (A. 81) has the following property: 



lim 

£->0 



1-e 



dy 



2e(y + l/3) , 20(-y + l/3) 



+ 



0. 



(A.82) 



i+, v(l-Z/)(3 + y) (l + y)(3-y) l-y\ 

Eq. (A.82) imphes that, although (in the absence of soft matrix element correc- 
tions [8]) HERWIG generates an incorrect angular distribution of soft gluons, the 
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total amount of radiated soft energy is identical to that computed in perturbative 
QCD. This is reassuring, and explains why HERWIG correctly describes physical ob- 
servables, eq. (A. 81) notwithstanding: any IR safe observable is not sensitive to the 
details of the angular distribution of soft gluons. 

Nevertheless, eq. (A. 82) prevents us from using dT,\Mc in the construction of the 
MC@NLO. In fact, dT,\MG a-cts as a local counterterm in eq. (4.22), and its soft 
behaviour (eq. (A. 81)) prevents the necessary cancellation. 

This problem can be solved in two different ways. The most straightforward 
procedure would be to define different showering variables, that would give short- 
distance cross sections with the correct soft limit. However, this would imply a 
modification of the existing MC code in one of its crucial parts, the showering al- 
gorithm. As already stressed, in this paper we prefer not to modify any substantial 
part of HERWIG. 

The second solution rehes on the observation that the angular distribution of 
soft gluons is irrelevant to MC@NLO physical results. We need a local counterterm 
to cancel the divergences of the real matrix element, but in the soft emission region 
the MC@NLO cross section is dominated by the pure MC contribution. On the other 
hand, for hard MC emission, we must have full control of the gluon emitted by the 
MC, since only in this way will the MC@NLO result match the NLO result. It follows 
that we can choose a form for dTi\yic identical to that we get by using eqs. (A. 72) 
and (A. 73) for non-soft emission, and a form that can act as a local counterterm for 
real matrix elements in the soft region. One possibihty is therefore the following 



d(f)3 



dcrqq{pi,P2) 



MC 



s \{l-x)^+ d{x,y) 



d(j)3 



MC 



s V(1-^)C- d{x,y) 



e{-y-Y,4x))g{x) 



{l-x){l + y) ^^ - m) ]Mmx^{^P,,xi'^P,) , (A.84) 



where ^ is a largely arbitrary function, such that 



lim g{x) = , lim g{x) = 1 . (A.85) 

>1 x—*xuz 



A simple choice is 



^ (l-n^))'" ^ / - ^Dz)/(1 - 5oz) X^z<X<l, 

r(^)2a + (l_J^(^))2a' |q X < X,, , 



(A.86) 



52 



where 



(A.87) 



The parameters in eqs. (A. 86) and (A.87) must be chosen in the following ranges: 
a > 1, < P < 1. The dependence of MC@NLO results upon these parameters will 
show the sensitivity of physical observables to the details of soft gluon emission. 

We point out that the procedure outlined above is not necessary when the matrix 
element is not singular in the soft limit. This happens for the partonic subprocesses 
qg W^W^q. In such a case, wc simply set ^ = 1 in eqs. (A. 83) and (A. 84). 

Finally, it is worth noting that the interesting features of eqs. (A. 81) and (A. 82) 
are of a general nature, since the hard production process is factorized into Adf^ 
in eq. (A. 81). Thus, the problem raised here is not specific to W'^'W" production. 
We expect that the angular distributions of soft gluons in HERWIG do not agree 
with those prescribed by the factorization theorems of QCD. Therefore, the local 
counterterm dT,\^c that we need in order to construct the MC@NLO will not coincide 
in general with the 0{as) HERWIG result. As in the present case, we expect that 
physical observables will not display any dependence upon the form chosen for c/SImc 
in the soft limit. In all cases, we must verify that the total amount of soft gluon 
energy given by the MC agrees with the QCD prediction, i.e., that the analogue of 
eq. (A.82) holds. 



B. Technicalities of MC@NLO 

In this appendix, we present formal proofs of some of the features of MC@NLO. 
We do so in the case of the toy model. However, we beheve that the corresponding 
properties hold in the case of QCD as well, as we shall point out explicitly; the toy 
model just allows us to simplify the notation. 

B.l Construction of Imc 

Our first aim is to give a precise definition of Imc{0,Xm{x)). We identify an MC 
event with the set of energies of the n photons emitted in the shower, {xk}k=i', 
the dependence upon x^ix) is understood. We define an MC@NLO event as the 
corresponding MC event plus the energy of the photon emitted at the NLO level: 
{x,Xk}k=i- In the case of soft emission, a; = 0, the MC@iNLO event simply coincides 
with the MC event. Any observable O will be in general a multi- valued function of 
MC@NLO events, 0{{x, Xk}l=^) = {d, . . . , Op}, with 1 < p < n + 1, the 0/s being 
values in the physical range of O. We then define 

p 

Af{{x, XkYU; O, (50) = ^ 9(0 + SO/2 - 0,)e(0,- - O + so/2), (B.l) 
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where O is any given value that the observable O can assume in its physical range. M 
is then just the number of entries in the bin (O — 50/2, O + SO/2) resulting from the 
MC@NLO event {x,Xk}k=i- Therefore 0<J\f<p<n + l. Starting from a given 
Xm{x), we generate N MC events, and we define 

_ 11^ _ 

lMc{0,x^ix)) = lim^hm-—J2^^^^^ 

~* 1=1 

rii being the number of photons, with energies x^'"*, emitted in event number i. Thus 
in general /mc(0, ^^(a;)) is the differential cross section that one obtains by running 
the MC starting at Xyi{x), defining O through all photons emitted in the shower, 
plus that emitted at the NLO level, normalized by the total number of MC events 
generated. 

If we assume that the MC is switched off, then from eq. (B.2) we obtain 

ly,c{0,x^{x)) = 6{0 - 0{x)), (B.3) 

0{x) being the single- valued function used at the NLO level to describe the ob- 
servable O. Eq. (B.3) can be proven by using eq. (B.2) in the computation of 

JdOlMc{0)f{0), for any arbitrary function f{0). With the same technique, we 
can also find a closed form for l-^ic{0,x^,i{x)) in the general case. We denote by 

) the probability that the toy MC emits n photons with energies 
xi > X2 > . . . > Xn- This probability can be obtained from eq. (3.14), and reads as 
follows 

Po = A(xo,Xm(x)), (B.4) 

Pn{Xi, ...,Xn) = a^A{xo, Xm{x))Q{Xm{x) - Xi) -Q{Xi - Xj+i), (B.5) 

i=l 

where in the case of n emissions we define Xn+i = xq- We also denote by the 
probability of having n emissions, regardless of the photon energies. Clearly, Pq = Pq- 
For n > 0, we can compute this quantity by integrating eq. (B.5), to obtain 

A{xo,x^{x)). (B.6) 

With eqs. (B.4) and (B.6) we can check directly that the toy MC conserves proba- 
bility: 

oo 

Po + ^Pn = l. (B.7) 

n=l 

We now introduce the functions On;j{x; Xi, . . . , x„) that return the multi- valued ob- 
servable O (for single- value observables, we simply have j = 1). Observing that 



P = 




^\dXi Pn{xi, 



: Xn) 



1=1 



a'' 
n! 



xo 



Q{z)_ 
z 
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the sum over all events in eq. (B.2) can be split into sums over events with a fixed 
number of emissions, we easily get the following: 

{0,x^{x))^Fo5{0-0{x)) 

oo „i / n \ p{x;xi,...,Xn) 

Xi , • • • , Xji ) S{0-0n;j{x;Xi,...,Xn)),{B.8) 

n=l / j=i 

where the integer p{x; Xi, . . . , x„) in general depends on the kinematics (this happens 
with jet observables, for example). Eq. (B.8) can be used, together with the explicit 
forms of the emission probabilities P„, to recover the results for /mc(2/) and /mc(^) 
given in eqs. (3.25), (3.26), (3.38), and (3.39); these equations have been obtained 
by direct computation, which constitutes a cross check of eq. (B.8). 

It should be noted that eq. (B.8) is obtained from eq. (B.2) without using the 
explicit form of the emission probabilities P„. Thus, both equations hold in the 
QCD case as well, with obvious formal substitutions (the MC@NLO event {x,Xky^^_^ 
gets replaced by an n + 1 kinematical configuration, and the information on parton 
flavours must be inserted). 

B.2 Total rate for exclusive observables 

We define an observable to be exclusive if its values can be obtained through a single- 
valued function, regardless of the number of the photons emitted in the shower. The 
hardest jet energy y, considered in section 3.4, is an observable of this kind, whereas 
the fully inclusive photon energy z is not. It is clear that in the exclusive case the 
function J\f defined in app. B.l has the following properties: 



Af{{x,xkrk=,;0,50)^ 0,1 
3!0 for SO^O such that : 



VO; 

A^({x,Xfe}Li;0,50) = i. 

Eq. (B.IO) holds for any given MC@NLO event. Then 



/ 



dOlMciO,XM{x)) = 1. 



(B.9) 
(B.IO) 



(B.ll) 



We can also get this results directly from eq. (B.8), using eq. (B.7), and the fact that 
the observable O is a single- valued function (and thusp(a;; xi, . . . , Xn) = !)• It is then 
apparent that eq. (B.ll) also holds in QCD. We can use eq. (B.ll) in eq. (3.20); we 
get 



/ 



dx 



+ 



r 

Jo 



dx 



dOIyic{0,Xyi{x)) 

dOI^c{0,l 
B + aV + a 



a[R{x) - BQ{x)] 



X 



B + aV + 
R{x) - B 



aB[Q{x) - 1] 



X 



X 



= CTtot- 



(B.12) 
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In the particular case of O = y, this result could be obtained from the results for 
da/dy given in sec. 3.4, by including the full xo dependence that has been neglected 
there. 

Since eq. (B.ll) holds in QCD, we conclude that eq. (B.12) also holds in QCD. 
This is in fact confirmed numerically by the results for the exclusive observables (such 
^ Pt^ ■> A^*"^^^) presented in sect. 4.7 (in the case that no cuts are apphed). 

B.3 Expansion to 0(as) 

We now consider the expansion to order a and of our formulae for the toy model 
and QCD MC@NLO, eqs. (3.20) and (4.22) respectively. We start with the toy model. 
Since we are not interested in C(a^) terms, only zero- and one-photon emissions from 
the MC are relevant. We have 



In eq. (B.13) we do not need to know the 0{a) term, since Imc{0,Xm{x)) is already 
multiplied by a in eq. (3.20). On the other hand 



The first term on the r.h.s. of eq. (B.14) is the contribution to the observable in the 
case that the MC does not emit any photons. The corresponding observable must 
be computed in the soft limit, with a coefficient which is the probability of no MC 
emission, that is, the Sudakov form factor for the range {xq, 1). The second term 
corresponds to the one-emission contribution. We stress that eqs. (B.13) and (B.14) 
could have been obtained directly from eq. (B.8). We can now substitute eqs. (B.13) 
and (B.15) into eq. (3.20). We get: 



7mc(0, Xm{x)) = 5{0 - 0{x)) + 0{a). 



(B.13) 




A(a;o, l)d{0 - 0{t)) + 0{a^) (B.14) 



(B.15) 
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a[R{x) - BQ{x)] 



dx S{0-0{x)) 



msub 



X 




(B.16) 
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We thus have 
da 
dO 



dx 



5{0-0{x))^^ + S{0 



X 



+ aB 



/ dx 

Jo 



Q{x) 



X 



0(0)) 

S{0 - O{0)) - 5{0 - 0{x)) 



B + aV-^ 

X 



+ 0(a'). (B.17) 



If we use eq. (B.17) to compute (O) = J dOO{da/dO), we get the NLO result of 
eq. (3.13), up to uncontrolled 0{a^) terms, and up to terms suppressed by powers of 
xq, which are numerically very small, and can be eliminated by letting xq — > 0. In 
the case of QCD, the role of Xq is played by some soft and coUinear cutoffs, which 
cannot be set to zero. Thus, in the real case we have power-suppressed corrections 
to the 0{as) result. Notice that MC's do not in any case control terms of this kind, 
and the dependence upon the cutoff is in general power-like. Since, as discussed in 
sect. 3.3.2, the MC@NLO and the MC behave the same with respect to the cutoff, 
the result in eq. (B.17) implies that the dependence upon the cutoff of the MC@NLO 
result beyond 0{as) is also power-like. 

Had we applied the same procedure using naive subtraction, eq. (3.19), we would 
have obtained the following (after letting — > 0): 



/ da 

\dd 



dx 



5{0 - 0{x))^^ + 5{0 - 0(0)) 



X 



+aB 



I 



B + aV- 



aB 



X 



X 



5{0 - 0{x)) - 5{0 - 0(0)) 



+ C»(a' 



(B.18) 



The second integral on the r.h.s. of this equation vanishes upon integration over the 
whole range of O, and the total rate (if meaningful) is then identical to that obtained 
from eq. (B.17). On the other hand, if we compute (O) the coefficient of 0{x) reads 
a{R{x) + BQ{x))/x, instead of the aR{x)/x that we get from eq. (B.17). This proves 
in general that the naive subtraction method suffers from double counting, as already 
shown in eq. (3.30) and (3.42). 

We now turn to the case of the QCD MC@NLO. The situation is completely 
analogous to that of the toy model. The NLO result can be recovered only if the 
0{as) terms of the expansion of /mc(0, 2), multiplied by the Born term, cancel the 
contributions of the terms dT,\Mc in eq. (4.22). The QCD case is however more 
involved, owing to the fact that, in the case of emission from initial-state partons, 
the Sudakov factor has a non-trivial dependence upon the parton densities. The 
analogue of eq. (B.15) in the QCD case reads as follows 



fMc(0,2) = 



L _ «s(/^Mc) f d^+ dz^ 
V 27r I 



27r 



27r 

d^^ dz_ 



■e(4-e+)Pi?(^+) 



1^+1 A*Mc) 



fiHl)( (MC) . 

Ja \Xi jfJ'Mc) 



AH2)( (Mr) 
Jc V 2 I 



■—■/'mc) 



/■(^^2)/ (MC) X 



+ 



5(0-0(2)) 
(B.19) 
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where the dots indicate two more terms of 0{as), identical to those appearing in 
eq. (B.19) except for the formal substitution 5{0 - 0(2)) S{0 - 0(3)). In 
eq. (B.19) we understand that the initial-state partons that initiate the shower have 
flavours a and b. The lower bound on the integration range, Qc, collectively denotes 
the cutoffs that prevent MC emissions from being soft and/or coUinear, and plays the 
same role as Xq in eq. (B.15). The arguments of the parton densities have already 
been discussed in app. A. 5. We also indicate explicitly the dependence upon the 
hard scale /Xmc, used by the MC in the shower; in principle, this scale depends on 
the kinematics of the emission. However, this is not important in what follows, and 
the notation we use simplifies the discussion. The Born term that gets multiplied by 
eq. (B.19) reads as follows: 

(B.20) 

The 0{as) terms in the product of eqs. (B.19) and (B.20) must cancel the contribu- 
tions of dTiablmc ill GQ- (4.22). It is apparent that a necessary condition for this to 
happen is xf'^'^'^ = xf''; this is precisely eq. (A. 60). However, the present proof adds 
one piece of information which could not be inferred from eq. (A. 60): in fact, that 
equation merely states the fact that the functional form of a;^''*^^ can be obtained 
through the procedure of event projection; however, it does not constrain the x\^'~''' 
variables to be equal to the x[^^ variables used in the event projection at the NLO 
level. This logical step is made only here. Once the condition on the Bjorken x's is 
fulfilled, we achieve the desired cancellation only if the short-distance cross section 
resulting from eqs. (B.19) and (B.20) match those in dTjablMc- The former have the 
form reported in eqs. (A. 72) and (A. 73), whereas we have chosen the latter as in 
eqs. (A. 83) and (A. 84); therefore, we seemingly have uncancelled 0{as) terms. How- 
ever, as argued in app. A. 5, the difference between the two forms cannot induce any 
effects on infrared-safe observables. This is in fact what we observe in practice (see 
sect. 4.7). In other words, we do not have event-by-event cancellation, but we do 
have cancellation on a statistical basis. We remind the reader that the uncancelled 
0{as) terms would disappear event by event if the MC had the correct angular dis- 
tribution for soft gluon emission. We point out that, owing to the presence of the 
cutoff in eq. (B.19), we still have a small uncancelled 0{as) contribution, which 
is the analogue of the last term on the r.h.s. of eq. (B.17). As in the case of the toy 
model, this contribution is suppressed by powers of the cutoff. Notice that, because 
of the simultaneous presence of soft and coUinear singularities in QCD, these powers 
can in general be multiplied by powers of the logarithm of the cutoff, but they remain 
negligible. 

We also mention the fact that the flavour structure of eq. (B.19) is non-trivial. 
For the cancellation at 0{as) to occur, we need to ensure that the flavour assignment 
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in the initial conditions for S events matches the flavour structure of the Born term, 
eq. (B.20). This is in fact guaranteed by the prescription described in sect. 4.5. 

We remark finally that the mismatch between /^-mc and /Xmconlo is formally of 
higher order. The two scales have fairly different physical meaning, and one cannot 
in general assume that /Imc — A*McaNLo- However, since the evolution of the parton 
densities in the MC is only leading-logarithmic, in the present paper we have made 
that assumption. A limited study of scale dependence has been presented in sect. 4.7. 

B.4 Resummation of large logarithms 

Here we give arguments to indicate that the MC@NLO resums large logarithms in 
the same way as the ordinary MC. As previoiisly, we start with the toy model, by 
reconsidering the results of section 3.4 for the variable z. We assume that eq. (3.37) 
is replaced by the more general form 

The function S is largely arbitrary, but it should at least have the following properties: 

• S has a power expansion in a, starting with an 0{a) term; 

• the leading behaviour at small z is a/z. 

The function in eq. (3.37) has these properties. Equation (B.21) is, however, more 
general; we can for example assume that S correctly takes into account next-to- 
leading logarithmic terms. Eqs. (3.38) and (3.39) now become 

I^ciz,!) = S{z,a), (B.22) 
Imc{z,Xm{x)) = S{z - x) + S{z,a)Q{xM{x) - z). (B.23) 



, R{z) - BQ{z) ^ R{x) - BQ{x) 

Gtot + a 777 ^ « / dx- 



We thus get 

Vc^^/^sub L zS{z,a) X 

Thus, since S{z, a) diverges as l/z for 2; — > 0, we have 

da 



(B.24) 



, ^atotS{z,a) + ... (B.25) 

where the neglected terms are either constant for 2; ^ 0, or vanish in this limit. 
Eq. (B.25) generalizes eq. (3.45). Formally, the 0{a) terms in atot produce a new 
tower of NLL terms when multiplied by the LL tower in S{z,a). But since Utot 
is factorized, this new tower of logarithms does not produce any change in the z 
spectrum of the emitted photons compared to the MC. This is in fact evident if we 
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normalize the z distribution to the total exclusive rate, which is <Jtot in the case of 
MC@NLO, and B in the case of MC: in both cases, we get S'(z,a). Furthermore, 
the terms neglected in eq. (3.45) do not lead to any logarithmic-enhanced behaviour. 
Thus, the leading and whatever subleading logarithmic behaviour in z of soft photon 
emission is the same in MC and MC@NLO evolution. 

Notice that eq. (B.25) imphes that the ratio of the number of soft photons to 
the number of showers is the same as in the ordinary MC. This fact, and the fact 
that the energy spectrum of the emitted photons is also the same, suggests that, for 
an arbitrary observable O that gets resummed through ordinary MC methods, the 
resummation performed by the MC@NLO has the same logarithmic accuracy as the 
MC. 

However, the preceding arguments are fairly closely linked to the details of the toy 
model, which makes it difficult to extend them to the QCD case. We therefore now 
give another argument, which can easily be extended to QCD; this argument works 
only for a given class of observables, which however includes most of the common 
measurable quantities. Namely, we consider the case of an exclusive observable O 
such that, when O ~ Os, any fixed-order cross section contains large logarithms 
log(0 — O^). Furthermore, we assume that, regardless of the number of photons 
emitted, at NLO and/or by the MC, when O ^ the energies of a// photons emitted 
arc constrained to be small; in other words, any hard photon emission moves O away 
from Os- This requirement is not very restrictive; still, there can be observables 
which do not meet it, and have large logarithms in the cross sections, but these 
logarithms will be multiplied by at least an extra factor of a relative to the case we 
treat. 

Our assumption implies that, when we consider Jmc(0,Xm(x)) for O — > O^, the 
energy x of the photon emitted at NLO is effectively a function of O: 



X 



= x{0) = k{0 -Os) + {{O - Osf) , (B.26) 



where we used the fact that x{Os) — 0. We can use eq. (B.26) to compute the NLO 
integrated cross section. We have 

JOs \"<^/nlo 

= / dx\ I dOS(0-0(x))^^^+ f dOS(0-Om(B + aV-—] 
Jo Uo. X Jo, \ X J ^ 



(B.27) 



where we have used eq. (3.13). With eq. (B.26) we obtain 



o^NLo(C') = B + aV+ [ dx 

Jo 



^e(.(o-o.)-x)-^ 

X X 



(B.28) 
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Consistently with eq. (3.4), we assume that 

R{x) = B + ^Rr,x'^ , (B.29) 

n=l 

and thus eq. (B.28) gives 

oo 

aNLo(O) =B + aB log {k{0 - O,)) + aV^ + a ^ — {k{0 - O,))" . (B.30) 

The first two terms on the right-hand side are the first two terms in the LL tower 
that is resummed by MC means. The third and fourth terms are NLL, the latter 
being also power suppressed. We point out that since the only large logarithm in the 
cross section has to be log(0 — O^), the coefficient k must be of order unity. 

We now consider the case of the MC@NLO. We denote the normalized pure MC 
result as follows: ^ 

a^c{0)= f dOI^a{0,l)- (B.31) 

We also notice that, owing to the fact that Xm{0) = 1, and assuming that Xm{x) has 
a power expansion in x around x — 0, we have (for O ~ Os) 

/ dOlMc{0,x^{x))= / (iO(/Mc(0,l) + NLL)e(«(0-0,)-x). (B.32) 

We can now repeat the same computation that led us to eq. (B.30), for the MC@NLO, 
eq. (3.20). In the case of § events, analogously to what happens for the contribution 
of the 0-body events in eq. (B.27), the restricted range in O does not constrain the 
integration range in x. Thus 



MO) = f 
Jo 



o. dO 



O-Mc(O) 



B + aV + aB l' dx^ 



/ 

Jo 



X 



(B.33) 



we therefore recover the pure MC result, up to 0{a) terms in the normalization. On 
the other hand, the restricted range in O does constrain the integration range in x 
in the case of H events, as shown in eq. (B.32). We obtain 



aaiO) = [ dO 

JOs 



dO 



= a (aMc(O) + NLL) {k{0 - Os)T , (B.34) 

n=l 

where Qn are the coefficients of the expansion of Q{x) around x — Therefore, 
eq. (B.34) does not contribute to the large-logarithm resummation for O — > Og. 
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Taking into account eq. (B.33), this allows us to conclude that the MC@NLO resums 
large logarithms in the same way as the ordinary MC does. 

It is clear that the argument given above works in the same way in the QCD case. 
One has to pay attention to the fact that, in QCD, soft and coUinear singularities 
exist. In general, the restricted range in O constrains the emissions to be either 
soft or coUinear, not necessarily soft and coUinear. This implies that in QCD each 
term of the power series that appears in eq. (B.34) can be multiplied by a power of 
log(0 — Os). However, this does not change our conclusion, that H events do not 
contribute large logarithms to the physical cross section. 



C. MC@NLO based upon slicing method 

In this appendix, we comment on the approach of refs. [14, 15], where a formalism is 
presented to construct an MC@NLO using the slicing method at the NLO level. This 
formalism is called $-space Veto in ref. [15]. It has the main objective and virtue 
of ensuring that only positive event weights are generated. However, we argue that 
this is at the expense of being able to recover the NLO results upon expansion in ots 
(this is because the resulting MC@NLO does not have an expansion in ctg), and a 
smooth matching between different kinematic regions. 

In the course of the discussion, we use the toy model, which is sufficient to high- 
light the main features of the method. We also propose a formula for an MC@NLO 
based upon the slicing method, which does not have any of the problems that affect 
the $-space Veto method, and fulfills the definition of MC@NLO given in sect. 2.1 
(but as a consequence is not guaranteed to generate positive event weights). 

The $-space Veto relies on the result for the NLO cross section obtained with 
the slicing method, eq. (3.9). The first step consists in fixing the parameter 5, by 
imposing the condition that eq. (3.9) does not contain any contribution proportional 
to 0(0). In other words, one seeks 5q such that^ 

S + a(Slog5o + y) = 0. (C.l) 

Thus 

B^aV ( B + aV \ 

With such a choice, eq. (3.9) becomes 

(0)_ = a + (C.3) 

JSo ^ 

Given the fact that the parameter 5o is not free any longer, being related to the 
Born and virtual contributions as shown in eq. (C.2), the uncontrolled 0{So) terms 

^This idea was first introduced, in a slightly different form, in ref. [12]. 
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in eq. (C.3) may be non-negligible. An improvement of eq. (C.3) is proposed in 
ref. [14], which uses a subtraction method in a restricted part of the phase space 
(this procedure is called "hybrid subtraction"): however, this improvement is of no 
importance in what follows, and will be neglected here (also, it is not used in current 
implementations of refs. [14, 15]). As anticipated, the reason for fixing the slicing 
parameter 5 can be seen in eq. (C.3); this equation does not contain any negative 
contributions, and this fact is used in refs. [14, 15] to define an MC@NLO with no 
negative-weight events. However, the resulting Sq is not analytic around a — 0, which 
seems to us to cause problems, as will be discussed below. 

The next step in the $-space Veto method is the introduction of a free parameter 
Sps, which fulfills the condition < ^ps < 1- In ref. [15], it is suggested that one set 
Sps equal to the boundary of the dead zone, in the case that HERWIG is used to 
construct the MC@NLO. In general, it seems difficult to prove that this condition is 
consistent with the requirement Sq < 5ps, especially since the boundary of the dead 
zone has a different (process dependent) shape from the region removed by slicing, 
and may depend on different variables. 

In the language of the present paper, the introduction of 5ps implies the definition 
of a modified slicing formula: 



(O) ,. = a dx 

\ / mslice I 







Oix)^^Qix - 5ps) + O(0)^e(a; - 5o)e(5ps - x) 

X X 



(C.4) 



Prom this equation, we can construct a toy MC@NLO as explained in sect. 3.3, 
paying attention to the fact that, in the $-space Veto method, the maximum energy 
available to the photon in the first branching is required to be 5ps (instead of 1) for 
class § events. Therefore 



,{0,Xy,{x))?^e{x-5,,) 

X 

+ /mc(0, 8p,)-^Q{x - 5o)e(5ps - x) 

X 



(C.5) 



Eq. (C.5) can now be used to compute the predictions of the $-space Veto method, 
in the same way that eq. (3.20) was used for our MC@NLO. We consider here the 
exclusive observable y, introduced in sect. 3.4. We can substitute eq. (3.26) directly 
into eq. (C.5); on the other hand, we need to use 

/Mc(y, <^Ps) = a^A(y, 5ps)e(5ps - y) (C.6) 

y 

instead of eq. (3.25). We find the following: 

( ^) = a^^Hv. S.s)ei6p, - y) dx ^ . (C.7) 

\dyJ^.nc. y J So ^ 
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We also find: 
o for y < Xe 



\dy J , 

\ C' / mslice 

o for y > Xg : 

dam 
dy 



',^A{y,xM) + a'^ r -R{x)A{y,xUx)) 



y Jsps 



X 



e(?/-5ps); 

(C.8) 



a 6(7/- 

y 



.5,,)+a'^e{x-'{y)-5,s) 

y 



dx 



<5ps 



X 



R{x)A{y,XM{x)) . 

(C.9) 

At ?/ ^ 0, we expect the result of the ordinary MC, eq. (3.24), to be reproduced, up 
to terms of higher logarithmic order, and up to multiplicative factors of higher order 
in a, i.e., we expect an equation similar to (3.35) to hold. Eq. (C.7) appears to fail to 
achieve this result, since it starts at (9(a^). However, 6q also depends upon a; thus, 
in order to reach firm conclusions, the integral in eq. (C.7) needs to be evaluated 
explicitly. To do that, we assume again that R{x) has a power series expansion in x, 
eq. (B.29). Eq. (C.7) then becomes 



dafi 



dy y 

' mslice y 

Using eq. (C.2) we now find 



'^^y^A{y,S,,)e{S,s-y) 



On — : n 



n=l 



(C.IO) 



das 
dy 



mslice 



a^A{y,6,sMS,s-y) 

y 



R. 



n / rn 



n=l 



n .'^PB-expl-n 



B + aV + aB log 5ps 
B + aV^ 



(C.ll) 



This equation has the same functional form in y as eq. (3.24), up to uncontrolled 
NLL terms (due to the fact that the second argument of the Sudakov form factor is 
5ps instead of 1), and formally starts at 0{a). However, it has the rather unpleasant 
feature that it cannot be expanded as a power series in a, since it is not analytic at 
a = 0. We therefore have the paradoxical situation that a quantity which is based 
on perturbation theory cannot be expanded in powers of the coupling constant. This 
is in fact due to the form of 5q, eq. (C.2). Still, if we assume that o^Sq is of O(a^), 
then eq. (C.ll) has the correct Sudakov behaviour for y ^ (this is also formally 
consistent with the necessity of neglecting powers of 5q in the slicing procedure; 
however, we believe that this is not justified here, since 5q is not a free parameter 
any longer). However, with the same assumption we conclude that the $-space Veto 
method has double counting problems, since for y < 5^^ the result of eq. (C.ll) does 
not coincide with the pure NLO result. On the other hand, for y > Sps the MC@NLO 
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and NLO results agree up to 0{a^) terms (see eqs. (C.8) and (C.9)). Therefore, the 
MC@NLO prediction has a discontinuity at y = 6ps- 

We now define an MC@NLO based upon the shcing method that does not have 
the problems outlined above. We can easily accomplish this task by noting that 
the slicing formulae can be obtained, at the NLO level, by making suitable approx- 
imations in the subtraction formulae. We start from eq. (3.9), which we re- write as 
follows 



Js 



0(x)^^ + 0(0) {B + aV + aBlogS) 

X 



(C.12) 



Prom now on, we do not write explicitly the terms of order 0{6) that are always 
neglected in shcing formulae. Notice that eqs. (3.9) and (C.12) in fact coincide up 
to 0{6) terms; the two equations would be exactly identical had we multiphed 0(0) 
by 1/(1 — 5). Using eq. (C.12) we now define an MC@NLO, following the procedure 
described in sect. 3.3. We obtain 



da 

mslice ^ o 



dx 



lyici0 1 2^m(2^)) 



a[R{x) - BQ{x)\ 

X 

Q{x) 



+/mc(0, l)[B + aV + aB\og5 + aB 



X 



(C.13) 



We stress that this formula has the same features as a "standard" slicing formula, 
such as eq. (C.12): the parameter S is free (and small), the integration over x does 
not extend down to zero, and the counterterm of the subtraction method is replaced 
by a term proportional to log 5. There is, however, one important difference. In 
eq. (C.12), the weights attached to 0{x) and 0(0) are divergent when 5 — > (and 
can now be negative). This does not happen in eq. (C.13) for the weights attached 
to Imc{0,Xm{x)) and /mc(0, 1). In the latter term, however, this property holds 
only if the integration over the phase space is performed before the limit 5 — > is 
carried out. This is inconvenient, because the unweighting procedure is performed 
in practice before the integration. In the spirit of the slicing method, we should 
integrate the term aBQ{x)/x in the weight attached to /mc(0, 1) analytically, which 
may or may not be possible, depending upon the form of Q{x). In order to simplify 
the discussion, we assume 

Q{x) = Q{x,^^,-x). (C.14) 
Using this form in eq. (C.13) we obtain 

aO J Jx X 

/ mslice ^ " 

+/mc(0, 1) [B^aV^ aB logx^^ad) ■ (C.15) 

In general, the weights attached to Iyic{0,XM_{x)) and /mc(0, 1) in this equation 
are not positive definite. However, it is clear that unweighting (in the sense of 
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generating events with weights of equal magnitude) can easily be achieved, since all 
the cancellations between large numbers occur under the integral sign. In this sense, 
the MC@NLO based upon the slicing method is similar to the MC@NLO based upon 
the subtraction method, eq. (3.20). However, while in eq. (3.20) the cancellation is 
achieved by numerical methods, the generalization of eq. (C.15) can be obtained from 
eq. (C.13) only through an analytical integration. This is similar to what happens at 
the NLO level; in the subtraction method, the residue of the soft/coUinear singularity 
is integrated numerically, whereas in the slicing method an analytical integration is 
performed. 

In the general case in which Q{x) does not give an easy analytical integral, we 
can use a mixture of analytical and numerical methods. Namely, we use the identity 



Q{x) = Q{xth - x) + Q{x) - Q{xtk - x) 



(C.16) 



in the last term of eq. (C.13), where Xth is an arbitrary constant. We readily obtain 

dx 



da 
dO 



X 



+Imc{0, l)iB + aV + aBlogXth + aB 



Q{x) - Q{xth-x) 



X 



. (C.17) 



We can now prove that eq. (C.13) has no double counting, and resums soft logarithms 
in the same way as standard MC does. We follow the procedures introduced in 
app. B.3 and B.4. Using eqs. (B.13) and (B.15), we get 



da 
dO 



— dx 

mslice *^ ^ 



6(0 - 0{x)) 



aR{x) 



X 



+ 6(0 - 0(0)) {B + aV + aB log d) 



+ aB / dx-^ ' 



X 



+ 0{a% 



d{0 - 0(0)) - 5{0 - 0{x)) e{x -5)- e{x - xo 

(C.18) 



which has to be compared to cq. (B.17). As in that case, wc recover the NLO result 
by considering JdO 0{da/dO), up to O(a^) and power-suppressed terms. Notice 
that we also have terms of 0{5) in eq. (C.18), consistently with the approximation 
upon which the slicing method is based. The distribution of the photon energies can 
be obtained as explained in app. B.4. We get 



mslice 



atot + a — r \-aB dx 

Zb[Z,a) Jmax(5,a:-i(z)) ^ 



, (C.19) 



which is the analogue of eq. (B.24). It is apparent that the same conclusions as in 
app. B.4 apply here for z — > 0. 
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Finally, notice that, by using eq. (B.ll), we obtain, in the case of an exclusive 
variable 

= atou (C.20) 

which is the result that we also find in app. B.2. We conclude that eq. (C.13) defines 
an MC@NLO that has the same formal properties as that defined in eq. (3.20). 



dx 



B + aV + aB\og5 + 



aR{x) 



X 
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